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Abstract 



An elementary introduction to the non-perturbative renormalization group is presented mainly in the context of 
statistical mechanics. No prior knowledge of field theory is necessary. The aim is this article is not to give an 
extensive overview of the subject but rather to insist on conceptual aspects and to explain in detail the main 
technical steps. It should be taken as an introduction to more advanced readings. 
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Chapter 1 

Wilson's renormalization group 



"What can I do, what can I write, 
Against the fall of night ?" 

A.E. Housman 

1.1 Introduction 

We give in these notes a short presentation of both the 
main ideas underlying Wilson's renormalization group 
(RG) and their concrete implementation under the form 
of what is now called the non-perturbative renormaliza- 
tion group (NPRG) or sometimes the functional renor- 
malization group (which can be perturbative) . Prior 
knowledge of perturbative field theory is not required 
for the understanding of the heart of the article. How- 
ever, some basic knowledge about phase transitions in 
Ising and 0{N) models is supposed. [H |2l [3] We shall 
mainly work in the framework of statistical field theory 
but when it will turn out to be illuminating we shall also 
use the language of particle physics. 

The beginning of this article will be rather elementary 
and known to most physicists working in the field of criti- 
cal phenomena. Nevertheless, both for completeness and 
to set up a language (actually a way of thinking at renor- 
malization group) it has appeared necessary to include 
it. The first section of this article deals with a com- 
parison between perturbative and Wilson's RG. Then, 
in the next section is presented the implementation of 
Kadanoff- Wilson's RG on the very peculiar case of the 
two-dimensional Ising model on the triangular lattice. 
This allows us to introduce the ideas of decimation and 
block spins and also those of RG flow of coupling con- 
stants and fixed point. The third section, which is also 



the heart of this article, deals with the "modern" imple- 
mentation of Wilson's ideas. The general framework as 
well as detailed calculations for the 0{N) models will be 
given. 

1.2 The perturbative method in 
field theory 

The idea behind perturbation theory is to consider an 
exactly soluble model, either the gaussian or the mean 
field "model", and to add, in a perturbation expansion, 
the term(s) present in the model under study and which 
are not taken into account in the exactly soluble model 
taken as a reference. [J UJ [S] For instance, in the "0''" 
model (which belongs to the same universality class as 
the Ising model) the partition function writes: 

Z[B] = J V(t>e-"'^'>'^+i ^'1' (1.1) 

with 

i/(0) = j d'x Q (V0)' + iro02 + 1^004^ (1.2) 

and B an external field and it is possible to take as a 
reference model the gaussian model: 

Zo = y"l?0e-^°('^'+^^'^ (1.3) 

where 

i/o(0) = / d'x (]- (V0)' -f ]-r^A . (1.4) 
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Z is then developed as a series in uo around Z^: 
Z = Jv^ (^1-^ ^\x,) 

(1.5) 

This expansion leads to the series of Feynman diagrams 
for the Green functions. 

The problem with this approach is that the "fluctua- 
tions" induced by the t/)^ term around the gaussian model 
are large. In the perturbation expansion they lead to in- 
tegrals — corresponding to the loops in the Feynman 
diagrams — of the form: 

/ d'^ql...d'^qL'[[{prop&g&tOY{q^)) (1.6) 

i 

where ^ 

propagator(gi) - . (1.7) 

(qi + Qr + ro 

These integrals are supposed to be cut-off at the upper 
bound by A which is an ultra-violet regulator. In the fol- 
lowing, it will be convenient to think at A as the (analog 
of the) inverse of a lattice spacing, lattice that would be 
used to regularize the field theory. In statistical mechan- 
ics, it is actually the other way around: the microscopic 
model is very often a lattice model whereas the field 
theory is an effective model only useful to describe the 
long-distance physics. 

If A was sent to infinity, the integral in Eq. (|1.6p would 
be generically divergent for d sufficiently large. This 
means that for A finite but large, the integrals are large 
and depend crucially on the value of A. This is very 
unpleasant for at least two reasons: 

(i) This invalidates the perturbation expansion even 
if uq is small. For instance, in the cj)'^ model 
and for the four-point connected correlation function 
G^\xi, . . . , X4) — {4>{xi) . . . 4>{x4))c, the one-loop ap- 
proximation writes in Fourier space at zero momentum: 

Gi'^ -^0 + (constant). ul I + . . . 

(1-8) 

This integral is divergent for d > 4 in the limit A ^ 00. 



(ii) The universal quantities (critical exponents, etc) 
are expected to be independent of the underlying lattice 
and thus, at least for these quantities, it is paradoxical 
that the lattice spacing {^-^ A~^) plays such a crucial role. 

Perturbative renormalization is the method that al- 
lows to reparametrize the perturbation expansion in such 
a way that the sensitive dependence on A has been elim- 
inated0 Then, the renormalization group allows to par- 
tially resum the perturbation expansions [6j and thus to 
compute universal behaviors. [H |2l El |5j 

Let us now make a list of questions that are not very 
often addressed in the litterature. Some answers are ex- 
plicitly given in this text. Some others require more 
thoughts. They are mainly there to nourish the reader's 
imagination... 

Ql: The occurence of ultra-violet divergences in field 
theory is often considered as a fundamental property of 
the theory. Thus, why do they play no role in the few 
known exact solutions of field theories or statistical mod- 
els ? For instance, in Onsager's solution of the two di- 
mensional Ising model no divergence occurs. This is also 
the case when Wilson's RG is implemented in field the- 
oretical models. 

Q2: Ultra-violet divergences are often said to be re- 
lated to the infinite number of degrees of freedom of a 
field theory (the value of the field at each point). But 
then, why does a classical field theory that also involves 
an infinite number of degrees of freedom show no diver- 
gence ? 

Q3: The answer to the last question often relies in 
the litterature on the fact that a statistical (or quantum) 
field theory involves fluctuations contrary to a classical 
field theory. Fluctuations are thus supposed to be re- 
sponsible for the divergences. The computation of the 
contributions of the fluctuations — the Feynman dia- 
grams — is thus often considered to be the reason why 
field theoretical techniques are relevant in statistical me- 
chanics. But there always exist thermal fluctuations in a 
statistical system whereas field theoretical techniques are 
most of the time useless in statistical mechanics. Then, 

^Let us emphasize that apart from the field renormal- 
ization, the whole renormalization process is nothing but a 
reparametrization. [B] 
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which types of fluctuations require field theory and which 
ones do not? 

Q4: Ultra-violet divergences are also often said to be 
related to the fact that we multiply fields at the same 
point (in a lagrangian) while fields are distributions the 
product of which is ill-defined. But what are the distri- 
butions in the case of the Ising model ? And since the 
interaction takes place between spins that are not on 
the same site but on two neighboring sites why should 
we take care about this difficulty ? 

Q5: In the cf)'^ theory for instance, the renormaliza- 
tion group flow of the coupling constant — given by the 
/3-function — is determined (in d = 4) by the UV di- 
vergences. But then why is the (IR stable) zero of the 
/3-function, that is the non-gaussian fixed point, useful 
to describe the infrared behavior of a field theory and in 
particular the critical behavior ? 

Q6: Why should we bother about the continuum limit 
in a statistical system — its ultra-violet behavior — for 
which on one hand there always exist a natural ultra- 
violet cut-off (such as a lattice spacing or a typical range 
of interaction) and for which on the other hand we are 
interested only in its long-distance physics ? 

Some of these questions will be answered directly in 
the following. For some others the reader will have to 
build his own answer. I hope that these notes will help 
him finding answers. 

1.3 Coarse-graining and effective 
theories 

There are two crucial remarks behind Wilson's 
method: [3H|9] 

(i) in general, we cannot compute exactly the con- 
tributions of the fiuctuations (otherwise we could solve 
exactly the model): approximations are necessary; 

(ii) the way fluctuations are summed over in pertur- 
bation theory is not appropriate since fluctuations of all 
wavelengths are treated on the same footing in Feynman 
diagrams. This is what produces integrals: at a given 
order of the perturbation expansion all fiuctuations are 
summed over[^ 

^ In quantum field theory, Feynman diagrams represent the 



Wilson's idea is to organize the summation over fluc- 
tuations in a better way. Note that because of remark (i) 
above, "better way" means "with respect to an approx- 
imation scheme"|j What is the idea behind Wilson's 
method of summation over the fluctuations? Before an- 
swering, let us notice that 

• in strongly correlated systems (e.g. close to a second 
order phase transition) the two relevant scales are (i) 
the microscopic scale a ^ — a lattice spacing, 
an intcrmolecular distance, the Planck length, etc 
— and (ii) the correlation length ^ (the mass(es) in 
the language of particle physics). These two scales 
are very different for T Tc and fluctuations ex- 
ist on all wavelengths between a and ^. In particle 
physics, ^ corresponds to the Compton wavelength 
of the particle {mc/H)^^ and a to the typical (in- 
verse) energy scale of the underlying "fundamental 
theory": lO""^^ GeV for a Grand Unified Theory or 
10^^ GeV for quantum gravity; 

• for the long distance physics and for universal quan- 
tities (magnetization, susceptibility, etc) the short 
distance "details" of the model have been completely 
washed out. This means that these "details" (exis- 
tence and shape of the lattice for instance) do mat- 
ter for the short distance physics but that they are 
averaged out at large distances: there must exist 
average processes that eliminate the microscopic de- 
tails as the scale at which we "observe" the system 
is enlarged 

Wilson's idea is therefore to build an effective theory 
for the long-distance degrees of freedom we are inter- 
ested in.[9l[8l[10l[TTl[12] This is achieved by integrating 
out the short distance ones. Since, at least for universal 

summation over probabiUty amphtudes corresponding to all possi- 
ble exchanges of virtual particles compatible with a given process 
at a given order. Note that these integrals are cut-oflF in the ultra- 
violet by A and in the infrared by the "mass" ro (see Eg. 111.61 1). In 
statistical mechanics, the mass is related to the correlation length 
5 by ro ~ (at the mean-field approximation). 

■^It is extremely rare that renormalization group enables to solve 
exactly a model that was not already solved by another and simpler 
method. 

*Note that this is true only for universal quantities. The criti- 
cal temperatures for instance, which are non-universal, depend on 
microscopic details such as the shape of the lattice. We shall come 
back on this notion in the following. 
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quantities, these short distance quantities do not matter 
crucially, it should be possible to devise approximations 
that preserve the physics at long distance. 

Actually, Wilson's idea is more general: it consists 
in saying that the "best" (approximate) way to study a 
subset of degrees of freedom of a system is to build an 
effective theory for them by integrating out the others. 
For instance, in molecular physics, one should build an 
effective hamiltonian for the valence electrons obtained 
by "integrating out" the core electrons (corresponding to 
high energy degrees of freedom). 

For the Ising model, this consists in integrating out 
in the partition function the "high energy modes" of the 
field (j){p) — those for which p G [A — dA, A] — and in 
computing the effective hamiltonian for the remaining 
modes. By iterating this procedure down to a scale fc, 
one should obtain an effective hamiltonian for the "low 
energy modes", those corresponding to p < k. The long 
distance physics, obtained for p ^ 0, should then be 
readable on the effective hamiltonian corresponding to 
A; ^ since no fluctuation would remain in this limit. 

Once again, let us emphasize that if we could perform 
exactly the integration on these "rapid" modes, we could 
iterate this integration and obtain the exact solution of 
the model. In most cases, this is impossible and the 
interest of this method, beyond its conceptual aspect, 
lies in the possibility to implement new approximation 
schemes better than the usual perturbation expansion!! 

Schematically, to implement Wilson's method, we di- 
vide (j){p) into two pieces: 4>^ (p) that involves the rapid 
modes p G [A/s, A] of (j){p) and (t>^{p) that involves the 
slow modes p € [0, A/s] : 

Z = I I?</)e-^[*'^'^l 

where K — {Ki ,K2,---) represents all possible coupling 
constants compatible with the symmetries of the sys- 
tem. Here, we have supposed that H involves all these 
couplings although the initial hamiltonian (that is, at 
scale A) can involve only a finite number of them. For 

^ Let us already mention that if Wilson's RG equations are 
truncated in a perturbation expansion, all the usual perturbative 
results are recovered as expected. 



instance, in the 0^ model all (initial) couplings Ki are 
vanishing, but those corresponding to the terms (V0)^, 
(jp and 0*. The integration of the rapid modes consists 
in integrating out the cj)^ field. In this integration, all 
the couplings Ki that were initially vanishing start to 
grow (this is why we have considered them from the be- 
ginning) but since we have considered the most general 
hamiltonian H (compatible with the symmetries of the 
problem), its functional form remains unchanged. Let 
us call K' the new coupling constants obtained after in- 
tegrating out (f)^. By definition of K': 

Z = J V(j)^e-"^^<-'^''^/'^ (1.10) 

with 

,K\A/s] ^ J ,K,A] (1.11) 

We thus build a series of coupling constants, each asso- 
ciated with a given scale: 

A -> K, 




K", etc. 

This method has several advantages compared with the 
usual, a la Feynman, approach: 

• There is no longer any summation over all length 
scales since the integration is performed on a mo- 
mentum shell, l^l S [A/s, A]. Thus, there can be no 
divergence and there is no need for any renormaliza- 
tion in the usual sense (subtraction of divergences). 

• The coupling constants Ki are naturally associated 
with a scale whereas this comes out in the perturba- 
tive scheme as a complicated by-product of regular- 
ization and renormalization. Wilson's method by- 
passes completely renormalization to directly deals 
with renormalization group. 

• The method is not linked with a particular expan- 
sion and there is therefore a hope to go beyond per- 
turbation expansion. 



6 



• The "flow" of coupling constants K ^ K' ^ A'" 
... is sufficient to obtain much information on the 
physics of the system under study. In particular, 
the notion of "fixed point" of this flow will play a 
particularly important role in statistical mechanics. 

1.4 Renormalization group trans- 
formations 

1.4.1 Blocks of spins 

As a pedagogical introduction, let us start by a sim- 
ple and illuminating example of Wilson's method im- 
plemented in x-space instead of momentum space and 
without having recourse to field theory. _1] [2] 





\ 1 / 
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\ / \ / 1\/ \ / \/ \ / 
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Figure 1.1: Partition of the triangular lattice in pla- 
quettes. The plaquettes are labelled by capital letters 
I,J,K,... and the spins inside the plaquettes are de- 
noted, in obvious notations, Sf, i = 1,2,3. The lattice 
of plaquettes is again triangular with a lattice spacing 
a\/3. 

We consider a triangular lattice with Ising spins to 
examplify block-spin transformations: 

H=-JY,S^SJ-BY,Si (1.13) 

<ij> i 

where Si are Ising spins: Si = ±1, the summation 
< ... > runs only on nearest neighbors and i? is a 
uniform magnetic field. The lattice is partitioned into 
triangular plaquettes labelled by capital letters I, J, ... . 
We call 5/, i = 1,2,3 the spin number i of the I-th pla- 
quette. As a first step, we separate the 8 configurations 



of the three spins S( of plaquette / into 2x4 configu- 
rations of (i) the block-spin Sj — ±1 — which is chosen 
here to be an Ising spin — and (ii) the four configura- 
tions, called cr"* , corresponding to a given value of Sj 
either -1-1 or —1. We choose (and this will be modified 
in the following) to define Sj by a majority rule: 

Si = sign(5( + S^ + Sl) = ±1 (1.14) 

Thus, for the four configurations of the S( compatible 
with 5/ = +1, we define the four variables a'^'^ by: 

cr]^ corresponds to ftt 

tjf"'" corresponds to tti 

L (1-15) 
a J corresponds to t|t 

(Tj^ corresponds to |tt 

For Si = —1, the cr"" correspond to the opposite con- 
figurations of the spins. Note that it will not be nec- 
essary in the following to compute the 5*/ in terms of 
the CT/^. The sum over all spin configurations in the 
partition function can be written as: 

The partition function can thus be rewritten as 

Z[B,T,n]^Y. ^ e-^I^-^*'^'^'"'"^!. (1.17) 
{5a{<±} 

where n is the number of lattice sites (n —^ oo corre- 
sponds to the thermodynamical limit) and a is the lattice 
spacing. We have also chosen to redefine the coupling 
constant J and the magnetic field B so that the pref- 
actor l/ksT of H in the Boltzmann weight is absorbed 
into the normalisation of these quantities. In this ex- 
ample, the lattice spacing of the lattice of plaquettes is 
aVS. The summation over the short distance degrees of 
freedom cr"^ can be formally performed 

Z[B, T,n]^Y. e-«'[^^'^'^'"/3--v^l (1.18) 

{Si} 

with, by deflnition of H': 

^-H'[Si,B,T,n/3,aV3] ^ ^-H[Si,af^ ,B,T,n,a] ^ -j^g-. 
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Let us make some remarks here. 

• We have chosen the majority rule to define the block 
spin Si so that it is again an Ising spin. The price to 
pay is that relation (|1.14|) between Si and the Sys 
is non linear. This is a source of many difficulties 
in more complicated systems than the Ising model 
on the triangular lattice. Moreover, it is difficult 
to generalize this rule to continuous, iV-component, 
spins. 

• The explicit computation of H', Eq. (|1.18p . shows 
that it involves infinitely many interaction terms, 
even if H involves only a nearest neighbor interac- 
tion. Thus, as it stands, the form of H is not stable 
under block-spin transformations. 

There is a solution to both problems. 

• We define Si through a linear transformation in- 
stead of the majority rule: 

(1.20) 

iei 

The Ising character is lost but the relation (|1.20p is 
simpler than (jl.l4p and can be generalized to other 
models. 

• We take for H a hamiltonian involving all possi- 
ble couplings among the Si compatible with the Z2 
symmetry {Si —Si for all i) of the Ising model 
(they will all be generated): 

<ij> <t.ij:> 

(1.21) 

+K3 SiS,SkSi + ... 

<ijkl> 

where <C ij ^ means summation over the next near- 
est neighbors. 

Now H ^ H{K, S„ n) where = (ifi, iCa, . . . ) repre- 
sents the set of all Z2-symmetric coupling constants as 
well as the magnetic field if necessary. For the initial 
hamiltonian: Ki^i = 0. Of course, this seems extremely 
complicated but it is the only possibility to have a form- 
invariant hamiltonian under block-spin transformations. 



The fact that all couplings are generated when fluctu- 
ations are integrated out simply means that even if an 
hamiltonian involves a finite number of couplings, all 
correlation functions, involving an arbitrary number of 
spins, are non-trivial. We shall come back on this point 
later. 

Let us finally remark that for n = cx), n/3 = 00, 
= 00, etc, and the "number of spins" remains iden- 
tical. However, the lattice spacing varies: a — > \/3a 
3a ~f . . . . Since, we shall look for fixed point hamilto- 
nian^ in the following, it will be necessary to rescale the 
lattice spacing by a factor l/VS after each summation 
over the rapid modes in such a way that we obtain, after 
summation over the CTj^'s, the same hamiltonian for the 
same system. We shall come back on this point later. 

We can now rewrite Eq. (jl.l9p for n — 00 as 

^-HIK\SimV3] ^ ^ g-H[A',5i,a°±,a]_ ^^ 22) 

This transformation, together with the rescaling of the 
lattice spacing, is called a renormalization group trans- 
formation0 Such a RG transformation 

• preserves the partition function Z and thus its sin- 
gularities and thus the critical behavior; more gener- 
ally, all thermodynamical quantities are preserved!! 

• maps a hamiltonian onto another hamiltonian (a 
system onto another system) in such a way that they 
have the same long distance physics; 

• consists in integrating out (averaging over) short 
distance degrees of freedom to obtain an effective 
hamiltonian for the long distance degrees of free- 
dom; 

^ Fixed point hamiltonian is meant in the usual sense. If 
H{K* , Si) is a fixed point hamiltonian this means that K* — > K* 
by summation over the it^^'s, Eg. HI. 191 1. The rescaling of the 
lattice spacing, which is equivalent to measuring all dimensionful 
quantities in terms of the running lattice spacing and not in terms 
of the (fixed) initial lattice spacing, is a necessary step to obtain 
fixed point hamiltonians. 

^We shall see in the following that this rescaling induces a 
rescaling of all coupling constants as well as of the magnitude of 
the spin, the so-called field renormalization. 

*If we wanted to compute correlation functions of the original 
spins, we would have first to couple the system to an arbitrary 
magnetic field (in order to be able to compute derivatives of Z 
with respect to the magnetic field Bi). This is a complicated task. 
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• can be summarized in a change of (infinitely many) 
coupling constants: K ^ K' . 

And now, two questions: 
Question 1: "Why is it interesting to integrate out the 
CTj^'s? Isn't it as complicated to integrate them out as 
would be the full integration over all degrees of freedom?" 

It is true that integrating out exactly the (T^^'s is of 
the same order of difficulty as calculating Z completely. 
However 

• the full calculation of Z contains much more infor- 
mations than what we want to obtain to get a satis- 
factory description of the critical physics. Moreover, 
for universal quantities, we guess that we shall be 
able to make rather drastic approximations as for 
the microscopic details of the model, that is the in- 
tegration of the short distance degrees of freedom, 
since they probably play a minor role; this opens 
the possibility of new approximation schemes; 

• the qualitative (or semi-quantitative) behavior of 
the RG flow of coupling constants K K' ^ 
K" . . . is enough to predict many non-trivial be- 
haviors occuring around a second order phase tran- 
sition. 

Question 2: "Why should we make a series of small 
block-spins (coarse-graining) instead of directly a large 
one?" 

This question, which is not independent of the first 
one, is a little subtle and requires some developments. 
Once again, if we were able to perform exactly the in- 
tegration over the (t"*'s, small or large blocks would 
make no difference. Thus, the problem comes from the 
approximations and is therefore not fully under control 
before precise calculations are performed. However, the 
general idea is not difficult to grasp. 

Let us call T{. ,p) the function that maps K = K'^^^ 
onto K'^p^ after p iterations of the RG transformations: 

K(p) =f {k^^\p^ (1.23) 

We, of course, have the property 

Kip) ^f(^K^''\p~r^ =f (f (i?(°),r) ,p-r) (1.24) 



and thus 

f (.,p) = f(f (.,r),p-r) . (1.25) 

This is called a self-similarity^ [13] propertyH If T was 
known exactly, this property would be trivially satisfied. 
However, once approximations are performed, it is gener- 
ically violated as is the case for instance in most pertur- 
bative expansions. 

Let us illustrate the concept of self-similarity on the 
simple example of differential equations.!^ We consider 
the trivial differential equation: 

y = -ey (1.26) 

with y{to) = yo- The solution is 

y = fit-to,yo)^yoe-'^'-'°K (1.27) 

Of course, / satisfies a self-similarity property which 
means that we can either (i) first integrate (|1.26p be- 
tween to and r to obtain y(r) = yr and then integrate 
again (|1.26p between t and t with y-j- as new initial con- 
dition or (ii) directly integrate (|1.26p between to and t: 

y{t) = fit - to, yo) = f{t - r, /(r - to, yo)) ■ (1.28) 

This is trivially satisfied by the exact solution (|1.27p 
since 

/(a, /3) = /?e-^" = /3e-^"e-'("-°) = /(a - a, /(a, /?)) 

(1-29) 

However, this property is violated at any finite order of 
the perturbation expansion in e of y{t). Let us show this 
at first order in e for which we, of course, obtain: 

y{t) = yo{l-e{t-to)) + 0{e') . (1.30) 

^ Something is said to be self-similar if it looks everywhere 
the same. In our case, the self-similar character comes from the 
fact that the functional form of the RG flow does not depend on 
the initial couplings K^'^^ since the same function T is used to 
transform i?*") into J?(P) or K^^^ into K^P\ This results in the 
fact that the right hand side of Eq. 11.251 1 is independent of r since 
the left hand side is. This independence is completely similar to 
the independence of the bare theory on the renormalization scale 
in perturbative renormalization (or of the renormalized theory on 
the bare scale). This is what allows to derive the Callan-Symanzik 
RG equations in the perturbative context. 
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This defines the approximation of order one of /: 

f^'\t-tQ,yo)^yoil-eit-to)). (1.31) 
We obtain at this order: 

/(I) (a - a, /(I) (a, /3)) = /3(1 - ea)(l - e(a - aj) 
= /(i)(a,/3) + e2/3a(a_a) 

(1.32) 

By comparing this resuh with Ea. (|1.29p . we find that 
self-similarity is satisfied at order e, as expected, but is 
violated at order e^. The problem is that this violation 
can be arbitrarily large if a (which represents t — to) is 
large. Thus, even if e is small, the self-similarity property 
is violated for large time intervals. This is true at any 
finite order of perturbation theory. This large violation 
comes ultimately from the fact that the perturbation ex- 
pansion is not an expansion in e but in e{t — tg). This 
is completely reminiscent of the perturbation expansion 
in field theory where the expansion is not performed in 
terms of uq but in terms of uq log A where uq is the bare 
coupling constant and A the cut-off (see Eq. (II. 8|) for the 
(j)'^ theory in d = 4). Reciprocally, it is clear that if 
a = t — to is small, so is the violation in Ea. p.32p since, 
in this case, both a and a — a are small. Thus, using 
perturbation expansions on small or, even better, on in- 
finitesimal time intervals preserves self-similarity at each 
step. In geometrical terms, this means that we can safely 
use perturbation theory to compute the envelope of the 
curve f{a,(3) — the field of tangent vectors given by the 
so-called /?- function — but not the curve itself. [6] The 
curve f{a,P) can only be obtained in a second step by 
integration of its envelope. 

The analogue for the RG is that small blocks will be 
under control. Coarse graining in this case respects self- 
similarity even when approximations are used while large 
ones lead inevitably to large errors as soon as approxi- 
mations are used. 

Before studying the structure of the RG flow, let us 
make two remarks about the RG transformations and 
their physical meaning. 

1.4.2 Two remarks concerning RG 
transformations 

The first remark is that it is still widely believed that 
the correlation length ^(r) is a measure of the typical 



size of clusters of spins having the same orientaion, that 
is of ordered domains (in the Ising case). As a conse- 
quence, it is believed that the divergence of ^ at Tj, is 
a consequence of the divergence of the size of these (so- 
called) naive clusters. The traditional metaphor is that 
at Tc there would exist oceans of up spins with conti- 
nents of down spins that would contain themselves lakes 
of up spins with islands of down spins, etc, with some 
kind of fractal geometry. This is wrong. It has been 
shown long ago that the distribution of cluster bound- 
aries does not scale at criticality. Rather, at a temper- 
ature Tp well below Tc the clusters of spins having the 
opposite sign of the spontaneous magnetization merge 
into a large percolating cluster. An important point is 
that, strictly speaking, no phase transition occurs at Tp 
since no local order parameter of the Ising model can be 
built out the spins in order to describe this transition: 
there is no singularity of the partition function at Tp. In 
fact, it is possible to construct clusters of spins that are 
critical at T^. These are the famous Fortuin and Kaste- 
leyn clusters. [2] They are used in the Swendsen-Wang 
algorithm of Monte Carlo simulations of the Ising model 
since they partially defeat critical slowing down. [15] 

The second remark is that the "microscope analogy" is 
often used to give an intuition of the physical meaning 
of the RG transformations. In this analogy, the coarse- 
graining implemented in the RG transformations would 
be similar to what occurs when the magnification of a 
microscope is varied. Let us imagine that we look at an 
image made of small pixels of definite colours (say blue, 
green or red). At a mesoscopic scale, the pixels are no 
longer seen and only a smearing of the colors of blocks 
of pixels can be observed. As a result, the "physics" 
observed would depend on the scale as in the RG trans- 
formations. This analogy has several virtues but also 
several drawbacks. Let us mention some. First our brain 
plays a crucial role for the color vision. From the three 
colors blue, green and red the cones in the retina are sen- 
sitive to, our brain is smart enough to reconstruct the 
impression of a continuous spectrum of colors. Although 
the analogy leads us to believe that our perception of 
the colors at a mesoscopic scale is a linear combination 
at this scale of the elementary colors of the pixels, this 
is not so. Second, in a RG transformation, there are two 
main steps (not to mention the final change of scale to 
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go back to the original lattice spacing). The first one is 
to build a stochastic variable for the block0 The sec- 
ond is to build an effective hamiltonian for this block 
variable by integration over short distance fluctuations. 
We can imagine that the first step is analogous to the 
superimposition of the colors of the different pixels. But 
then what is the analog of the second step? The laws of 
classical electrodynamics for the propagation of light do 
not change from one scale to the other. Let us repeat 
here that the effective hamiltonians for the block vari- 
ables in the Ising model are extremely complicated: they 
involve all powers of the fields and not only interactions 
among nearest neighbors There is no analog for this 
step in the microscope analogy although it is the crucial 
one from the RG point of view. In fact, things go almost 
the other way around. Whereas the electromagnetic field 
emitted by several pixels is the linear superposition of the 
field produced by each of them, the /3-function in quan- 
tum field theory that gives the evolution of the coupling 
constant with the scale is a measure of the deviation to 
the trivial rescaling invariance (in the case of quantum 
electrodynamics). Thus, although the microscope anal- 
ogy can be useful it should be employed with some care 
(and a grain of salt). 

Let us now show how linear RG transformations can 
be implemented. This will allow us to prove a simple 
relation about the behavior under RG transformations 
of the two-point correlation function (SiSj) and thus on 
the correlation length. [U [9] 

1.4.3 Linear RG transformations and 
behavior of the correlation length 

Instead of the majority rule, we consider a linear trans- 
formation between the spins of a plaquette and the block- 
spin. [T] The simplest idea is to take a spatial average (not 
a thermodynamic one) 

S'j^Ws! (1-33) 

^"This can be performed either by a majority rule as in Eo. 111. 1411 
or by a linear relation as in Eg. 111.341 1. 

^^Once the continuum limit has been taken and continuous RG 
transformations are implemented this means that the effective 
hamiltonians involve all powers of the field and of its derivatives. 



where s is the "linear size" of the block, that is s is the 
number of spins per block. In our example of the trian- 
gular lattice, d = 2 and s = ^/S. As we already said, we 
shall also need to perform a rescaling of all dimensionful 
quantities (in order to find fixed points). Thus we take: 

5, = ^^5/ (L34) 

where A(s) is a function that will be determined in such 
a way that we find a fixed point. This relation among the 
stochastic variables S( and Sj leads to relations among 
their thermodynamic averages. The most important one 
is the two-point correlation function: 

{SiSj) =1^5,5,e-^[^^^-l 

{Si} 

g2d E G'^^H^i^Xj) 

(L35) 

where, by definition 

G(2)(x„f,) = G(2)(ry,i?) = {S,S,) (1.36) 

is the two-point correlation function of the spins St. If 
the correlation length is large compared with the size of 
the plaquettes and if we consider two plaquettes / and 
J such that their distance is very large compared to a: 
Tij — \xit^i — Xj^j\ 3> a, then G'^'^^ {xi,Xj) does not vary 
much for i £ / and j £ J. Thus, in this case: 

G^^Hx,,x,)~s^^G^^Hx,,x,). (1.37) 

Therefore, close to the critical temperature and for dis- 
tant plaquettes: 

~A2(s)(5/5/). (1.38) 

The important point is that (SiSj) is also the two-point 
correlation function of a Z2-invariant magnetic system. 
The only difference with (SiSj) is that it is computed 
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with the set of coupUngs K' instead of K. We thus 
obtain: 

(5/5j> = G(2)(r,j,i?') (1.39) 

where G'^' is the same function as in Eq. (|1.36p and r/j 
is the distance between the plaquettes / and j[3 We 
thus deduce: 

G(2) {rij , i?') - (s) G(2) [n, ,K) (1 .40) 

for sufficiently distant plaquettes. 

Let us now explain the precise meaning of this rela- 
tion. Let us suppose that we are given a new model on 
a triangular lattice with a set of couplings K' . In prin- 
ciple, we can compute the correlation function G^^-* of 
the spins of this new system. Our claim is that this cor- 
relation function is identical to the correlation function 
of the block-spins Sj of the original system. Of course, 
to compare the two functions we have to say how to 
compare the distances r^- between spins in the two lat- 
tices. Our calculation shows that what we have to do 
is to measure all distances in the length unit intrinsic 
to the system, that is in units of the lattice spacing of 
each system. This means that the quantities and rjj 
that appear in Eq. (jl.40p are pure numbers that must be 
numerically diSerent 

rij = rij/s (1.41) 

since they correspond to the same "distance" but mea- 
sured in two different units: a for the original system 
and a' = sa = a\/3 for the coarse-grained system. It is 
important to understand that measured in an extrinsic 
length unit, like meters, these distances are indeed the 
same: r/ja' — (rij / s).{sa) — rija whereas they are "dif- 
ferent" — they correspond to different numbers — when 
they are measured in the length unit intrinsic to each sys- 
tem. Put it differently, the dimensionful distances r/ja' 
and Tij-a, measured in a common length unit, are equal 
whereas the dimensionless distances r/j and r^, mea- 
sured in terms of the lattice spacing of each system, are 

Let us point out here a subtlety. This statement is not fully 
rigorous since the original spins are Ising spins whereas the block 
spins Sj are not. The correlation functions (5^ S^) and (SjSj) are 
therefore not computed exactly in the same way since the summa- 
tion over the configurations of Si and of Si do not run on the same 
values. In fact, after several blocking iterations, the spins that are 
summed over become almost continuous variables and the afore- 
mentioned difficulty disappears. 



different. We shall put a bar or a tilda on dimensionless 
quantities to distinguish them. 

Let us emphasize that the value in an extrinsic unit 
like meters of, say, the correlation length is almost mean- 
ingless. From a physical point of view, the only relevant 
measure of the correlation length is in units of the lat- 
tice spacing. The difhculty in our case is two- fold. First, 
to write down a field theory, it is necessary to perform 
the continuum limit a — > 0. It is therefore convenient 
to rescale the position vectors by a factor a before per- 
forming this limit (thus, as usual, [x] =length). This is 
consistent with the fact that the vectors Xi of the lattice 
have integer components that label the position of the 
sites of the lattice and are therefore also dimensionless 
whereas, in field theory, the modulus of x represents the 
distance to the originli^l Second, since we have to con- 
sider several systems with different lattice spacings, it 
will be convenient to work in the continuum with lengths 
measured in units of the running lattice spacings. 

Let us therefore define dimensionless quantities as 



where r and ^ are the dimensionless quantities that will 
be convenient once the continuum limit will be taken, 
that is in the field theory formalism. Eq. ()1.40|) that in- 
volves only dimensionless quantities can then be rewrit- 
ten: 

G(2) c^X^{s)G^^\f,K). (1.43) 

Let us give a concrete example of the meaning of this 
relation. In three dimensions and at large distances: 
rij ^ 1, a typical form of the two-point correlation func- 
tion is: 

(5,5,)=G(2)(r-0~^ (1.44) 

with ^ the correlation length in units of the lattice spac- 
ing a and f = fij. We can use the same formula as in 

^^It will also be convenient to rescale the spin-field by the 
appropriate power of a: Si — > 0(5;) so that the gradient 
term (VcA)^ comes in the hamiltonian of the field theory with 
a dimensionless pre- factor (chosen to be 1/2 for convenience): 
H = J d'^x (i{V(/))2 -I- U{<f>)). We find from this equation that 

_ — 2 _ d-2 

l<f>(x)] = [x 2 ] SO that the rescaling involves a factor a 2 . 
Note that the original variables Si are dimensionless since Si = ±1. 
The function G^^^ in Eq. I|1.40| I is therefore also dimensionless. 
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Eq. (|1.44p for the correlation function of the block-spin 
system: 

-f'/f' 

{SiSj) ^ ^—^ . (1.45) 
r ' 

We thus obtain 

_f/£ 

G^'^ir^^')-"-^^' (1-46) 
and, by comparing with Eq. p.40p . we find that 

= - and A(s) = _ 47) 
s 

By comparing all these relations we find that 

• the dimensionful correlation lengths of the original 
and of the block-spin systems are identical: it is a 
RG-invariant. Thus, the dimensionless correlation 
lengths decrease as the coarse-graining scale s in- 
creases. This means that the coarse-grained systems 
are less correlated than the initial one and that the 
correlation length decreases linearly with the scale. 
Since the correlation length behaves as a power law 
close to the critical temperature 

C^(r-r,)-^ (1.48) 

parametrizing G*-^-* in terms of the (reduced) tem- 
perature 



or in terms of the correlation length is "equivalent". 
Thus, saying that f decreases with the block size s is 
equivalent to saying that the running reduced tem- 
perature t{s) increases with s: the coarse-grained 
system is "less critical" than the original one. We 
call relevant a parameter that increases with the 
scale s. 

• the reduced temperature t is one particular coupling 
among all the couplings K. We shall explain in the 
following why the form of G*^^^ given in Eq. (|1.44p 
and in which only the correlation length appears is 
valid at large distance that is why all other couplings 
in K play no role at large distances. 



• if we combine two RG transformations of scale si 
and S2 we must obtain the same result as a unique 
transformation of scale S1S2 (this is self-similarity). 
This clearly implies that 

A(si)A(s2) = A(siS2) . (1.50) 

It is straightforward to show that the only solution 
of this equation is a power law. The example above 
shows that the exponent of this power law is directly 
related to the power law behavior of G'-^-' (f) at T = 
r„ Eq.([L47|. 

1.5 Properties of the RG flow: 
fixed points, critical surface, 
relevant directions 

The RG flow takes place in the space of hamiltonians, 
that is in the space of coupling constants K. We now 
study this flow. One of its nice properties is that, with- 
out specifying any particular statistical system, very gen- 
eral informations on second order transitions can be ob- 
tained from it by only making very natural assumptions. 

At Tc, ^ ~ 00 (^ = 00) at a second order phase transi- 
tion. Thus the point J?'^"^ is mapped onto K^^^^ under a 
RG transformation for which ^' = 00 again (the block- 
spin system is also critical). We define the critical surface 
as the set of points K in the coupling constant space for 
which ^ = 00. For a second order phase transition only 
one parameter needs to be fine-tuned to make the sys- 
tem critical (the temperature for instance). Thus, the 
critical surface is of co-dimension one. It is stable under 
RG transformations. 

If we now consider a system described by a point K^^"^ 
such that ^ < 00 then ^' = ^/s and the block-spin sys- 
tem, being "less critical", is described by K^^'' which is 
"further away" from the critical surface than If 
we iterate the blocking process, we obtain points K'-^\ 
K^^\. . . that will be further and further away from the 
critical surface. 

We shall consider in the following the continuum limit 
of the Ising model and this will allow us to perform con- 
tinuous RG transformations. We call in this case the 
set of points Kg, s £ M., a RG trajectory. To different 
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Figure 1.2: Schematic representation of the RG flow in 
the space of coupHngs, i^i, . . . . This space is infinite 
dimensional and the critical surface, defined by the set 
of points for which the correlation length is infinite is 
of co-dimension one. Under RG transformations K 
K' K" . . . , the critical surface is stable: if K = Kc 
is on the critical surface, then K' , K", etc are also on the 
critical surface. A point away from the critical surface is 
mapped onto another one "further away" from it. 

Kg on the same RG trajectory correspond systems that 
are microscopically different (this means at the scale of 
their own lattice spacing) but that lead to the same long- 
distance physics since they all have the same partition 
function. Let us now make the fundamental hypothesis 
that must be checked on each example: |9] 

Hypothesis: For points in a (finite or infinite) domain 
on the critical surface, the RG flow converges to a fixed 
point K*: K* =f (^i?*,s). 

This domain is called the basin of attraction of the 
fixed point K* . Under this hypothesis, the typical topol- 
ogy of the flow on the critical surface is summarized in 
Fig ll.3l We have called "physical line" in this figure the 
line on which the temperature alone is varied. It is not a 
RG trajectory. Reciprocally, a RG trajectory does not, 
in general, correspond to any transformation doable on 
a physical system by a human being. It is only a map- 
ping that preserves the partition function without any 
connection to a physical transformation. 

All systems that belong to the basin of attraction 
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Figure 1.3: Schematic representation of continuous RG 
trajectories on the critical surface. The flow converges 
to the fixed point K* . For a given model, the "physical 
line" corresponds to a change of the temperature. It is 
not a RG trajectory. For a more precise description of 
the RG flow, see FigEH 

of K* belong to the same universality class since they 
all have the same long-distance physics. Note that K* 
depends on the choice of RG transformations T{.,s). 
Apart from being fixed for this particular choice of RG 
transformations, this point is nothing else than a partic- 
ular critical point. 

1.5.1 Scaling relations — linearization of 
the flow around the flxed point 

The existence of an attractive (in the critical surface) 
fixed point is sufficient to explain universality since, inde- 
pendently of the starting point Kc on the critical surface, 
all RG trajectories end at the same point K* . However, 
universality holds also for systems that are not right at 
Tc (which is anyway impossible to reach experimentally) 
but close to Tc- It is thus natural to assume that the 
flow is continuous in the vicinity of K* . In this case, 
starting at a point K close to the point and on the 
same physical line, the RG trajectory emanating from 
K remains close to the one emanating from Kc during 
many RG steps before it diverges from the critical sur- 
face, see Fig ll.4l It is easy to estimate the typical value 
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Figure 1.4: Schematic representation of two continuous 
RG trajectories corresponding to the same model for two 
different temperatures. The trajectory starting at Kc is 
on the critical surface (T = Tc) and the other, starting at 
ifo, is slightly away from it. There exists a RG trajectory 
emanating from the fixed point K* and which is not 
on the critical surface. It is an eigendirection of the 
RG flow corresponding to the relevant direction. Note 
that we could parametrize the coupling constant space 
in such a way that the axis denoted represents the 
temperature. It is not necessary that this axis coincides 
with the relevant eigendirection of the RG flow at K* 
but it is necessary that it has a non vanishing projection 
onto this axis since the temperature is for sure a relevant 
parameter. Note that this is non trivial in an infinite 
dimensional space. 

of s for which the RG trajectory diverges from the crit- 
ical surface. 

As long as the running (dimensionless) correlation 
length ^(s) — ^/s remains large, the system behaves as 
if it was critical and the representative point K(s) must 
be close to the critical surface. When the running cor- 
relation length becomes of order 1, the coarse grained 
system is no longer critical and K{s) must be at a dis- 
tance of order 1 of the critical surface. More precisely, at 
the beginning of the flow, K{s) moves towards K* (by 
continuity). It remains close to it as long as ^(s) remains 
large so that the memory of the initial point K is largely 
lost. Finally, it departs from K* when s ~ ^. Another 



and more precise way to state the same result is to say 
that the running reduced temperature t(s) is of order 1 
when s ~ ^: 

t{s) - 1 for s - f . (1.51) 

We shall see in the following that the hypothesis of the 
existence of a fixed point together with this relation 
are sufficient to predict the existence of power law be- 
haviors for many thcrmodynamical quantities with crit- 
ical exponents that are universal and that satisfy "scal- 
ing relations" independently of any specific microscopic 
modeJ. [9] Clearly, to obtain these relations, only the 
vicinity of K* is important since after a few RG steps, 
all RG trajectories emanating from points close to criti- 
cality are in the vicinity of this point, see Fig |1.4l This 
will allow us to linearize the RG flow around K* . 

For the sake of simplicity, we assume in the following 
that s can take continuous values (we work in the con- 
tinuum, that is the continuous limit, a — > 0, has been 
taken). We also suppose that the RG transformation 
with s = 1 is the identity (no block of spins) and that 
the composition law is such that a RG transformation of 
parameter si followed by another one of parameter S2 is 
equivalent to a transformation of parameter siS2- 

As already emphasized, to obtain approximations of 
the RG flow that are under control, it is preferable to per- 
form a series of infinitesimal RG transformations rather 
than directly a transformation with s large. It is thus 
necessary to study the differential form of these trans- 
formations. There are two ways to do so: 

(i) The first way consists in comparing two RG trans- 
formations of parameters s and s + e performed from an 
initial set of couplings Kin. 

Ks+e -Ks=f (ifi„, s + e) - f (Xin, s) (1-52) 

and thus 



This formula is not convenient since it is expressed in 
terms of the initial couplings Ki^ and, as a consequence, 
in terms of a parameter s that can be large (and thus 
for RG transformations that will not be under control 
once approximations will be implemented). Using ap- 
proximate expressions of T together with this expression 
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would lead to violations of the composition law of RG 
transformations. 

(ii) The second way consists in comparing two sets of 
running couplings differing by an infinitesimal RG trans- 
formation: Ks and i^T^^^i+g) and in expressing the result 
as a function of : 



T(K,A 



T(K,A 



and thus 



ds 



dT 
97 



(1.54) 



(1.55) 



In this second formulation, the evolution at one step — 
that is between s and s + ds — depends only on the 
couplings at this step — that is at s — and not on those 
at the previous steps: the evolution of the couplings is 
local in the space of couplings. 

It is important to notice here several things. 

• The composition law of RG transformations will be 
automatically satisfied if K{s) is obtained by inte- 
gration of this expression, even if an approximate 
expression of T is used. This comes from the fact 
that, by construction, K{s') is computed in terms 
of K(s) by composing the infinitesimal group law 
between s and s': this is the meaning of integrating 
Eg. ([135)1 . 



• Eq. (|1.55p leads naturally to the logarithmic evolu- 
tion of the couplings with s, contrarily to Eq. (|1.53p . 



In Eq. (|1.53|) . the variations of T had to be known 
for only one set of couplings i?in but for all s. This 
was the real problem of this approach since we are 
interested in large values of s. In Eq. ()1.55|) only the 
variation of T for s = 1 have to be known. The price 
to pay is that this must be known for all values of 
Ks- In perturbation theory this is not problematic 
as long as the running coupling used in the pertur- 
bation expansion remains small all along the how. 
Of course, if Kg converges to K* such that the cou- 
pling(s) of the perturbation expansion at the fixed 
point is large, perturbation theory runs into trouble. 



We define 



(1.56) 



which gives the evolution of the couplings of the model 
with the scale. Note that (3 (^Ks^ does not depend ex- 
plicitly on s since the right hand side of Ea. (|1.56p is 
evaluated at s = 1 and thus does not depend on this 
variable^ The only dependence on s of this function is 
implicit and comes through the dependence of the cou- 
plings Kg on s. By definition of the fixed point K* 



= . 



(1.57) 



We thus obtain 

dK, dK* 



ds 



ds 



(3 (Kg) - (3 [k*) 



Cli\ o 



where, by definition 



5Kg ^Kg-K* 



and 



is a matrix 



dA 
AK 



s,3 



K' 



Thus, in the neighborhood of K" 



M5Kg 



(1.59) 



(1.60) 



(1.61) 



M is not symmetric in general. We suppose in the fol- 
lowing that it can be diagonalized and that its eigenval- 
ues are reaf^ (this must be checked on each example). 
We, moreover suppose that the set of eigenvectors {ci} 
is a complete basis: 



and 



Mei = Xi 



5Kg 



with 



vAs) e. 



(1.62) 



(1.63) 



^*In perturbation theory, the fact that the /3 function of the 
marginal coupling is cut-off independent and thus scale indepen- 
dent is a consequence of the perturbative renormalizability of the 
model. 

^^It can happen that some eigenvalues are complex. In this 
case, the RG flow around the fixed point is spiral-like (focus) in 
the corresponding eigendirections. 
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Under these hypotheses, we obtain: 



S ; = } ^ Vi[S)M Cj 



ds 



s a 



and thus 

dvi{s) 



Xi Viis) 



(1.64) 



(1.65) 



We conclude that around the fixed point, the RG flow 
behaves as power laws along its eigendirections (but if 
Ai = in which case it is logarithmic). There are three 
possibilities: 

• Ai > and Vi{s) y when s y which means that 
the flow in the direction goes away from K* . et is 
called a relevant direction and Vi a relevant coupling. 
As we have already seen, the reduced temperature 
is relevant since the system is less and less critical 
along the RG flow which means that t{s) increases 
with s. 

• Ai < and Vi[s) \ when s /* which means that the 
flow in the direction approaches K* . Ci is called 
an irrelevant direction and Vi an irrelevant coupling. 

• Ai and Vi is said to be marginal. It is neces- 
sary to go beyond the linear approximation to know 
whether it is relevant or irrelevant. The flow in this 
direction is slow: it is logarithmic instead of being a 
power law. This is important since it is the case of 
the renormalizable couplings in the critical dimen- 
sion (that is d = 4 for the Ising model). 

Physically, we expect a small number of relevant direc- 
tions since, clearly, there are as many such directions as 
there are parameters to be fine-tuned to be on the criti- 
cal surface (the co-dimension of this surface) . For second 
order phase transitions there is one coupling to be fine- 
tuned to make the system critical: the temperature (for 
instance). There is actually one more parameter which 
is relevant but which is also Z2 non-invariant: the mag- 
netic field. All the other directions of the RG fiow are 
supposed to be irrelevant (this can be checked explicitly 
once a specific model is given). 



In the following, we shall use the language of mag- 
netic systems although our discussion will be completely 
general. We shall suppose that together with the mag- 
netic field, there is only one other relevant direction. We 
shall show that this implies scaling laws for all thermo- 
dynamical quantities with only two independent critical 
exponents. 

1.5.2 The correlation length and the 
spin-spin correlation function 

Let us start by studying the critical physics of a model 
at zero external magnetic field having only one relevant 
coupling. We call vi this coupling and Ai the eigen- 
value of the flow at K* associated with vi. We order 
the other eigen-couplings V2jV3,. . . in such a way that 
> A2 > A3 . . . . Physically, the reduced temperature is 
expected to be a relevant parameter. It does not nec- 
essarily correspond to the relevant eigendirection of the 
RG fiow but its projection onto this eigendirection is non 
vanishing. For the sake of simplicity we ignore this sub- 
lety that does not play any role in the following and we 
assume that vi = t. We have seen, Eq. (|1.40|) . that for 
large r: 

G(2) (f, k) ~ X-'is) G(2) (J, i?,) with X{s) = s'l^ . 

(1.66) 

From these relations we can deduce two others, one at 
Tc and another one away from Tc- 

• For i = and in the neighborhood of K* : 

G(2) (f , 0, i;2, 1^3, . . . ) - s-'g(2) 0, S^^«2, S^^«3, . . . ) . 

(1.67) 

Let us now suppose that we integrate out all fiuctuations 
between scales o and r. This amounts to taking s — f: 

G(2) (f , 0, z;2, . . . ) =^ f-'^Gf") (1, 0, f^^«2, . . . ) ■ (1-68) 

Since A2 < 0, we obtain that for r suflaciently large: 

f^'^V2 <C 1 and thus 

G(2)(f,0,i;2,...)-^"'G(^' (1,0,0,...) . (1.69) 



From the definition of the anomalous dimension 

1 



G(2)(r) 



oc 



(1.70) 
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we find that 

9 = d-2 + r]. (1.71) 

Note that the d— 2 part of 6 is purely dimensional: it cor- 
responds to the engineering dimension of the spin field. 
The rj part corresponds to a dynamical contribution. It 
can be proven rigorously and for many field theories (as 
the (j)'^ theory) that f] > 0. This means that the fluc- 
tuations contribute to decrease the correlations of the 
system with the distance. 

• For t 0, ^ < oo and we can therefore integrate 
all fluctuations between scales a and ^ by taking s — ^. 
As explained above, Ea. (|1.5ip . the running temperature 
at this scale must be of order 1 since the coarse-grained 
system at this scale is no longer critical. We thus deduce 
for large correlation lengths 

t{s = 0-t^^' ^1 . (1.72) 

From the definition of the critical exponent v: 

^r^t-" (1.73) 

we find that 

^=^. (1.74) 

The behavior of the correlation function close to the crit- 
ical temperature follows from Eas. (|1.66ll.72p and from 
?»1: 

G(2) {f,t,v,,...) ^ r'G(2) (^:,^p^«2p^...().75) 
^r^G(2) (J, 1,0...) 

- r-'f (0 . (1.76) 

We can see on this relation that close to Tc and at large 
distance the correlation function is no longer a function 
of infinitely many coupling constants but of only one 
parameter, the correlation length. One can also see that 
the Yukawa-like form of the correlation function that we 
have considered in Eq. ()1.44p has the right form. 



1.5.3 Scaling of the correlation function 
— Relations among exponents 

We now couple the system to a uniform magnetic field. 
In a RG transformation: 

BEs^ = Bj:^fi-B.j:S,. (1.77) 

Thus 

Bs^ s'^-"''^ B . (1.78) 

We call 

Xs=d-6/2='^±^ (1.79) 

the magnetic eigenvalue. Since 77 is always smaller than 
d + 2 (even for d = 1), > so that the magnetic field 
is a relevant variable. 

Relation (jl.66p can be trivially generalized in the pres- 
ence of a magnetic field and to any correlation function. 
We now consider the magnetization per spin which is the 
one-point function: G*^^-* — m = (S). Clearly, it behaves 
as: 

m{t, B,...)c^ s-^^^mis^^ i, B, . . . ) . (1.80) 

Once again, we can obtain several relations among ex- 
ponents by considering the scaling of physical quantities 
at and away from Tc- 

• For i = and B 7^ 

to(0, B, . . . ) ~ s"'*/2m(0, s^^ S, . . . ) . (1.81) 

and, by taking s such that s"^^ B ~ 1, we obtain 

m(0,B,...) ~ B''/2^^m(0,l,...) 

cx S^/* (1.82) 

by definition of the exponent 6. We thus find: 

(1—2 + 77 

• For t <0 and i? = O we find: 

m{t, O, . . . ) ~ s~^/^m{s^' i, O, . . . ) . (1.84) 
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and by taking s — ^: 

m(t,0,...) ~ ^-^/2m(l,0...) 

oc (T, - Tf (1.85) 

by definition of the exponent j3. Tlius we find 

d- 2 + 77 



(3 = v- 



(1.86) 



• For the susceptibility x = dm/dB, at t 7^ 0, we 
obtain: 

x{t,B,...) s-'/^-^m{s>^- t,s^-B,...) 

~ s^B-e/2^(s^ii,s^BB,...) . (1.87) 
By taking B = and s = ^, we find: 

x(i,o,...) ^ e^--«/2x(i,o...) 



Thus, by definition of 7 

7 = 1/(2 - 7?) . 



(1- 



(1.89) 



Finally, for the exponent a of the specific heat, the 
calculation is subtler and requires to consider the free 
energy. One finds: 



= 2-vd . 



(1.90) 



To conclude, we have found that the hypothesis of the 
existence of a fixed point in the RG fiow is sufficient to 
explain: 

• universality since the critical exponents depend only 
on the RG fiow around the fixed point and not on 
the point Kc representing the system when it is crit- 
ical; 

• the scaling behaviors of the thermodynamical quan- 
tities such as the magnetization, the susceptibility, 
the correlation fimction, etc; 

• the relations existing between critical exponents; 

• the irrelevance of infinitely many couplings and the 
fact that the scaling of the correlation length with 
the temperature drives the scaling of many other 
thermodynamical quantities. 



Two remarks are in order here. First, for second order 
transitions, only two exponents are independent, 1/ and 77 
for instance. Second, universality is a much more general 
concept than what we have seen here on critical expo- 
nents. It is possible to show in particular that the RG 
theory enables to understand why it is possible to keep 
track of only a small number of coupling constants in 
field theory while these theories involves infinitely many 
degrees of freedom and thus, a priori, infinitely many 
couplings. 

1.5.4 The example of the two- 
dimensional Ising model on the 
triangular lattice 

This model is very famous as an example where the 
block-spin method a la Kadanoff can be implemented 
rather easily. [16j We shall not repeat the explicit calcula- 
tion of the RG fiow that can be found in most textbooks 
on this subject. [T] We shall only give the main ideas that 
will be relevant for our purpose. 

As already explained above, the idea is to partition the 
lattice in triangular plaquettes, to build an Ising block- 
spin Si for each plaquette by a majority rule, Eq. (|1.14p . 
and to integrate out the fiuctuations inside the plaque- 
ttes compatible with a given value of Si. We have al- 
ready said that the implementation of this idea imposes 
to take into account all possible Z2-invariant couplings 
among the spins. From a practical point of view it is 
thus necessary to perform truncations of this infinite di- 
mensional space of coupling constants. The approxima- 
tions usually performed consists in keeping only some 
couplings and projecting the RG fiow onto this restricted 
space of couplings. We shall see that almost the same 
idea is used in the field theoretical implementation of the 
NPRG. The simplest such truncation consists in keep- 
ing only the nearest neighbor interaction, that is Ki , in 
hamiltonian (|1.2ip as well as the magnetic field B. For 
concreteness, let us give here the result of a RG transfor- 
mation with s — \/3 on the triangular lattice for a small 
magnetic field when only Ki and B are retained: 



Bs ^ W 



3e-ifi + e3^i 



(1.91) 
(1.92) 
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A fixed point is trivially found 



Kl = 0.336 (1.93) 
B* = . (1.94) 

Tfie exponent v can be computed as well as S and the 
result is: 

1/ = 1.118 , (5 = 2.17. (1.95) 

All the other exponents can be found using the scaling 
relations among them. The exact values, given by On- 
sager's solution are 

iy=l , (5=15. (1.96) 

Let us notice that these RG results can be systematically 
improved by keeping more and more couplings. They 
drastically and rapidly improve with the next orders of 
approximations. [T^ 
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Chapter 2 

The non-perturbative renormalization 
group 



"O insense qui croit que je ne suis pas toi !" 

V. Hugo 

2.1 Introduction 

All the different implementations of the non- 
perturbative renormalization group (NPRG) rely 
on Kadanoff- Wilson's ideas of block spins, coarse 
-graining and effective long-distance ttieories. However, 
they can substantially differ as for the way they are 
implemented. In the framework of field theory, there 
exists two main formulations : the Wilson (also called 
Wilson-Polchinski) approach^ HI [HI [HI [Tl] and the 
"effective average action" approach^ EQI Ell [221 EH 
[24, 25, 26, 27, 28,|29l[30l[31, 32, 33j. We shall deal with 
the second one which is not the best known, probably 
for historical reasons. Since it is nevertheless interesting 
to have an idea of the Wilson-Polchinski formulation, 
we start by this approach although we shall not study 
it in detail. 

2.1.1 The Wilson-Polchinski approach 

We shall work in the context of statistical field theory (at 
equilibrium). This means that we shall not deal with a 
minkowskian metric (this brings its own difhculties) and 
that we suppose a continuum description of the systems 
under study. The microscopic physics is supposed to cor- 
respond to a scale A in momentum space which is — up 



to a factor of unity — the inverse of a microscopic length 
(a lattice spacing, an intermolecular distance, etc). The 
partition function is thus given by a functional integral : 



Z[B] 



dlJ-cM) exp 



t/(0) 



(2.1) 



where cL^Ca is a (functional) gaussian measure with a 
cut-off at scale A: 



d^CA = 'D4>{x) exp 



1 



with (in momentum space) 

and C is the usual free propagator: 

Cip) = ■ 
p~ + r 



(2.2) 
(2.3) 

(2.4) 



The "cut-off" function 0^ is a step function in p-space 
starting at A and smoothened around A on an interval 
of typical width e, see Fig |2.1l If e = 0, the quadratic 
part of the hamiltonian becomes the usual gradient and 
"mass" term, cut-off at scale A: 



(2.5) 



The role of e is to smoothcn the sharp cut-off at A which 
is conceptually simple but technically unpleasant. 
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i-e,(p) 
1 



A-e 



- P 



Figure 2.1: A typical cut-off function in the Wilson- 
Polchinski approach. 

We want to implement the block-spin idea in our field 
theory framework, that is to separate (f)p = 4>{p) into 
"rapid" and "slow" modes (with respect to a scale k). 
The slow modes will play the role of block-spins while 
the rapid ones will correspond to fluctuations inside the 
blocks. It is convenient to work in Fourier space where 
the derivative operators are diagonalized (and the cut-off 
simple) . 

We define 

0p = 0p,< + 0p,> (2.6) 

and associate 



Ckip) 



(2.7) 



is the sum of 0^ and 



^p,, ^ Ca{p) ^ Ck{p) 

It is important to notice that 
(j)p^^ for allp: it does not coincide in general with (f>p,^ on 
[0, k] and with (pp^-^ on [fc, A]. The meaning of 0^ and 
comes from their propagator Ck{p) and Ca{p) — Ck{p) 
respectively. A beautiful identity allows us to rewrite the 
original partition function in terms of (j)^ and 0^ , (p) 
and C\{p) — Ck{p) and to perform (at least formally) the 
integration on the rapid modes: 



df^CA (0) = dfiCk (</-< ) dficA-c, (0> ) 
Let us show it on a one dimensional integral: 



/ = 



+ 00 



(2.8) 



(2.9) 



where x is the analogue of and 7 of Ca- We now define 



y + z 
a + P 



(2.10) 



y and z are the analogues of 4. 
a and f3 the analogues of Ca 



and 4>^ respectively and 
Ck and of Ck- Then 



Jcx 



+ C30 



dy dz e 



(2.11) 



We prove this identity in the Appendix, section [3721 Ac- 
tually, we are not only interested in gaussian integrals 
and our result can be trivially generalized: 

+00 /'+00 



(2.12) 

This result can be generalized straightforwardly to func- 
tional integals since it is a property of the gaussian inte- 
gralsQ It becomes 

dMCA(0)e-^^(«(x y'dMc.(0<) dA^CA-c.(0>). 

.e-/v'(0<+0>) (2.13) 

Thus, by performing formally the integration on (p^ , we 
define a running "potential" Vk at scale k: 

^-lv.{4>J^ /d^CA-c.(0>)e-^^(^<+^>) (2.14) 



with, by definition, Va — V, the initial potential. This 
definition leads to: 



dficMe-^'''^'^<^ ■ (2.15) 



Let us emphasize that Vk is called a potential because 
we do not have included in it the quadratic derivative 
term. However, generically, as soon as fc < A, Vk in- 
volves derivative terms with, moreover, any power of the 
derivatives of (p^ ■ 

It is possible to write down a differential equation for 
the evolution of Vk with fc: this is the Wilson-Polchinski 
equation derived in the Appendix section [3. 3.1 1 [9l [8l [12] 
It is a possible starting point for a non-perturbative for- 
mulation of the RG. The one we shall use is mathemat- 
ically equivalent to this one although more convenient 
in many respects when approximations are used. Be- 
fore going to this second formulation, let us make some 
remarks here. 



^We first consider an A^-dimensional gaussian integral and then 
take the limit — > 00. 
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(j>^ at scale k represents approximately a spatial av- 
erage of (j) over a volume of order k~'^. It is not 
a thermal average. For fc = 0, 0p=o,< represents 
only what is (improperly) called the magnetization 
mode: / d'^x 4>{x) which is not the magnetization. It 
has a highly non trivial probability distribution and 
the true magnetization is the thermal average of this 
mode. Thus (j)^ at scale k ^ is not a precursor of 
the order parameter: it is still a stochastic variable 
whose physical interpretation is not so trivial. 

The flow of "potentials" Vk{4>^) does not contain all 
the informations on the initial theory : the corre- 
lation functions of the rapid modes cannot be com- 
puted from this flow. It is necessary to first couple 
the system to an arbitrary "source" B{x) and to fol- 
low the flow of this term to reconstruct the corre- 
lations of the initial fields in the whole momentum 
range p G [0,A]. Thus the information is splitted 
into two different kinds of terms: on one hand, the 
fc-dependent hamiltonians which give rise to a flow 
for all the couplings involved in these hamiltonians 
and, on the other hand, the source term. Fortu- 
nately, much informations about the theory (e.g. 
critical behavior) can be obtained from the flow of 
hamiltonians alone. It is nevertheless an open ques- 
tion to know if the difficulties encountered within 
Wilson-Polchinski's method are related to the fact 
that the informations on the Green functions are 
not contained in the flow of hamiltonians. 

The effective hamiltonians Vk (and the correspond- 
ing Boltzmann weights) are highly abstract objects ! 
One should remember in particular that the RG 
transformations do not correspond to any transfor- 
mations that a human being can perform on the 
system. This is a purely theoretical idea that more- 
over will be useful only when approximations are 
used. 

The flow equation on Vk was written more than 
thirty years ago but was not much used in actual 
physical problems before the mid 90's (for excep- 
tions see, Golner, Newman, Riedel, Bagnuls, Bervil- 
lier, Zumbach, etc). [HI [Ml IMl ISi ISZl [351 ISSl SO] 
There are three main reasons for this strange mat- 
ter of fact. First, perturbation theory was ex- 



tremely successful for 0{N) models (as well as in 
particle physics) and the need for the NPRG was 
not obvious in many situations. Second, (renor- 
malized) perturbation theory was believed to cor- 
respond to a controlled approximation whereas ap- 
proximations performed in the NPRG framework 
seemed uncontrolled. However, this is completely 
wrong. Renormalized perturbation series are not 
convergent. [m [5] They are asymptotic series, at 
bestH For the 0{N)' model at Z — 1 loops, the 
coefficient in front of of the 4-point correlation 
function behaves at large I as II {—aY with a a real 
number. Thus, even in cases where many orders 
of the perturbation expansion have been computed, 
they are useless as such and resummation methods 
of the renormalized series are required (Pade-Borel, 
conformal-Borel, etc). In many cases, these resum- 
mation techniques fail to produce converged results. 
As for the NPRG, there is no general theorem about 
the convergence of the series of approximations that 
are used. However, from the few results yet ob- 
tained, it seems that this method has good conver- 
gence properties. [32[|42ll43l [4ll45l [46| It is however 
too early to draw any firm conclusion on this ques- 
tion. Third, it seemed that the anomalous dimen- 
sion was crucially depending on the choice of cut-off 
function 9^ that separates the rapid and the slow 
modes whereas it should be independent of it. This 
was especially important when studying the 0{N) 
models in two dimensions where it seemed impossi- 
ble to reproduce the perturbative results obtained 
from the non linear sigma model. This difficulty is 
very simply overcome in the "effective average ac- 
tion" approach. [32 11471133] 

Let us finally mention two other "psychological" diffi- 
cuhies related to NPRG. 

• The NPRG equation on Vk can be truncated in per- 
turbation theory. This, of course, enables to recover 

^The 0{N) models are completely exceptional in this respect 
since they are the only ones for which it has been proven that the 
series of the /3-function is Borel-summable in ci = 3 (in the so-called 
zero momentum massive scheme). In all other cases, either this is 
not known or it is known that the series are not Borel-summable. 
For QED, this is not yet a problem because the smallness of the 
coupling constant, the fine structure constant, ensures up to now 
an apparent convergence of the perturbative results. 
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the usual perturbation expansion (what else could 
it lead to ?). However, the way it was implemented 
most of the time in the 70's did not allow to retrieve 
the two-loop resultsll It is still (erroneously) widely 
believed for this reason, even by "specialists", that 
Wilson's method does not work at two-loop order 
and beyond! 

• Vk{(t>) involves infinitely many couplings contrarily 
to perturbation theory that involves only the renor- 
malizable ones. For this reason, it is widely believed 
that the recourse to numerical methods is unavoid- 
able in the NPRG approach whereas it is not in 
the perturbative one. This is not fully correct for 
two reasons. First, even in perturbation theory the 
RG flow cannot be integrated analytically in gen- 
eral. Second, even in the NPRG approach, very 
crude approximations, involving only very few cou- 
plings, often lead both to analytically tractable com- 
putations and to highly non-trivial non-perturbative 
results. [32] 

Let us now turn to the other implementation of the 
NPRG formalism. 

2.1.2 The effective average action 
method 

Many formal results about the NPRG method as well 
as some "physical" results have been obtained within the 
Wilson-Polchinski approach^ However, the revival of 
Wilson's ideas as well as their concrete implementation 
in the last fifteen years is largely linked with the devel- 
opment of an alternative, although formally equivalent, 
formulation first promoted by C. Wetterich at the end 
of the 80's under the name of effective average action 
method.m IMISSllMllllllHllinilSniinilSl] in prac- 
tice, this has allowed to compute in a reasonable way 
the anomalous dimension 77 and, more importantly, to 

''These calculations did not correspond to a series expansion of 
the exact NPRG equation on Vj.. They enabled to retrieve the 
one-loop results easily but became very cumbersome beyond one- 
loop. 

^See the impressive and inspiring worlds of Bagnuls and Bervil- 
lier about the formal aspects of Wilson's RG, as well as their crit- 
icisms of the perturbative approach. 1481 1491 1491 15U| 



study the physics of the 0{N) models and of many oth- 
ers, in all dimensions, including two.[32l [33] Moreover, 
the whole scheme is more intuitive, allows to retrieve 
very easily the one-loop results both in 4 — e and 2 + e 
dimensions and in the large N limit. This has convinced 
many specialists of the subject to work with this formal- 
ism. 

Block-spins, coarse-graining, Legendre trans- 
form, etc 

The original Kadanoff- Wilson's idea is to perform a 
coarse-graining and thus to map hamiltonians onto other 
hamiltonians at larger scales. The hamiltonians thus ob- 
tained are the hamiltonians of the modes that have not 
yet been integrated out in the partition function, that is 
4>^ . As already emphasized, these hamiltonians are very 
abstract objects^ Instead of computing this sequence 
of hamiltonians, we can compute the Gibbs free energy 
r[M] of the rapid modes (that is (j>^ ) that have already 
been integrated outlf] Thus, the idea is to build a one- 
parameter family of models, indexed by a scale k such 
that [5^ 

• when fc = A, that is when no fluctuation has been 
integrated out, the Gibbs free energy rfe[M] is equal 
to the microscopic hamiltonian0 

Tk=A[M]^ H[q}^ M] . (2.16) 

• when fc = 0, that is when all fluctuations arc in- 
tegrated out, Tk=o is nothing but the Gibbs free 
energy of the original model : 

Tk^o[M]^T[M]. (2.17) 

*It is impossible to get easily any physical information from it 
except at "mean field-like" level: a functional integral has still to 
be performed. 

^The Helmoltz free energy is _F = — fegTlog Z[B]. It is a func- 
tional of the source B(x). The Gibbs free energy is obtained by 
a Legendre transform from F and is a functional of the magne- 
tization M{x). It is the generating functional of the one-particle 
irreducible (IPI) correlation functions. 

^Let us emphasize that at the mean-field approximation, the 
Gibbs free energy of the system is identical to the hamiltonian. 
Ea. ll2.16l l is an exact version of this statement (remember that no 
fluctuation is taken into account at the mean- field level). 
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Thus, as k decreases more and more fluctuations are in- 
tegrated out. The magnetization at scale k is therefore a 
precursor of the true magnetization (obtained at fc = 0) 
and the free energy Tk, also called the effective average 
action, a precursor of the true free energy F (also called 
the effective action). 

Let us notice two points. First, k plays the role of an 
ultra-violet cut-off for the slow modes 0^ in the Wilson- 
Polchinski formulation (analogous to A in the original 
model). It plays the role of an infrared cut-off in the 
effective average action method since Fj. is the free en- 
ergy of the rapid modes. Second, the slow modes play 
a fundamental role in the Wilson-Polchinski approach 
whereas they are absent of the effective average action 
method which involves only the (free energy of the) rapid 
modes. As a by-product, we shall see that all the in- 
formations on the model (RG flow, existence of a fixed 
point, computations of correlation functions, etc) are 
contained in a single object : Ffc[M]. This is a rather 
important advantage of this method compared to the 
Wilson-Polchinski one. 

Now the question is to build explicitly this one- 
parameter family of Ffc. The idea is to decouple the slow 
modes of the model in the partition function. A very 
convenient implementation of this idea is to give them 
a large mass.[32| In the language of particle physics, a 
large mass corresponds to a small Compton wavelength 
(h/mc) and thus to a small range of distances where 
quantum fluctuations are important. A very heavy par- 
ticle decouples from the low energy (compared to its 
mass) physics since it can play a role at energies below 
its mass threshold only through virtual processes. These 
processes are themselves suppressed by inverse powers of 
the mass of the heavy particle coming from its propaga- 
tor. In the language of critical phenomena, the "mass" 
term rcfp' /2 in the hamiltonian corresponds to the devi- 
ation from the critical temperature: r (x T — Tc (at the 
mean field level, at least). Thus a large "mass" r corre- 
sponds to a theory which is far from criticality (^ ~ a), 
that is where thermal fluctuations are small. 

Therefore, the idea is to build a one-parameter family 
of models for which a "momentum-dependent mass term" 



has been added to the original hamiltonian : 'WI] 

Zk[B] = I Vcl>{x)exp (^-H[^] - AHk[cf>] + J B<j)^ 

(2.18) 

with 

AHk[(t>]^^J^Rk{q)^g4'-g- (2.19) 

The function Rk{q) — called the cut-off function from 
now on — must be chosen in such a way that 

• when fc = , Rk^oil) ~ identically (Vg) so that : 

Zk^o[B] = Z[B] . (2.20) 

This will ensure that the original model is recov- 
ered when all fluctuations are integrated out, see 
Eq.dllTl). 

• when fc = A, all fluctuations are frozen. This will 
ensure that relation (j2.16|) is satisfied. Giving an 
infinite mass to all modes q G [0, A] freezes their 
propagation completely and we must therefore take 

Rk=Aiq) = oo yq. (2.21) 

An approximate, but convenient way to achieve this 
goal is to choose a function Rk=A not infinite but of 
the order of A^ for all momenta. 

• when < k < A, the rapid modes (those for which 
l^l > k) must be almost unaffected by Rk{q) which 
must therefore almost vanish for these modes : 

Rk{\q\>k)~{) (2.22) 

On the contrary, the slow modes must have a mass 
that almost decouple them from the long distance 
physics. 

Remembering that Rk is homogeneous to a mass square, 
it is not difficult to guess its generic shape, at least if we 
require that it is an analytic function of q^ , see Fig l2.2l 
We shall discuss in the following some convenient 
choices for Rk but let us first define precisely what Fj, is. 
In principle, having defined by Eq. (|2.18[) . the Leg- 
endre transform of log Zk should be unambiguous and 
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Figure 2.2: A typical cut-off fimction in tlie effective 
average action approach. 



sliould lead to F^. Let us follow this program and see 
that there is a subtlety. We thus define 



Wk[B]= log Zk[B] 



(2.23) 



which is the Helmoltz free energy, up to the —ksT term 
that plays no role in what follows. The Legendre trans- 
form of Wk is defined by 



F',[M] + Tyfe[B] = / BM 



(2.24) 



where the magnetization M{x) is, by definition the av- 
erage of 0(a;) and is therefore : 



M{x) 



6B{x) 



(2.25) 



Of course, for fc ^ 0, Rk 0, Wk W and thus 
FJ,, F= Gibbs free energy of the system. However, 
it is easy to show that F^[A/] ^ H[M] contrarily to 
what is expected, see Eq. (|2.16|) . This comes from the 
AHk=A term, which is large. Thus, we prefer to work 
with a modified free energy where the Rk term has been 
subtracted in Ffc.lSj We define 



Tk[M] + Wk[B] = / BM 



Rk{q)MgM^g . (2.26) 



The Rk term in Eq. (|2.26p does not spoil the limit fc ^ 
since in this limit it vanishes Vg. Let us now show that 
Ea. (|2.26p is the correct definition of F^ leading to the 
limit Tk=A[M] = H{M), Ea. (l2T7)l . 



2.1.3 An integral representation of 
and the limit A; ^ A 



We start from the definition of Zfc, Eq . (|2.18p and from 
the definition of Ffc, Eq. (|2.26p . We deduce by differenti- 
ation (see. 13.3.21 ) 



B^ — 



5Tk 
5M^ 



+ / Rk{x-y)My 



(2.27) 



Thus, by substituting Eq. ((2?26)l and Eq. ((2?27|) into the 
defintion of Wk we obtain : 



-rfc[j\/] 



exp 



Vcfyexp 



5Tk 



((/.,- M,) Rk{x - y) {^y - My)^ (2.28) 



If we choose a function Rk {q) that diverges for all q as 
fc — !■ A then, in this limit : 

exp (^-IJ i<Px - M^) Rk{x - y) {(f)y - My)^ - 5{(I)~M) 

(2.29) 

that is, it behaves as a functional Dirac delta. Therefore, 



Vk[M]^ H[(t) = M] as fc^A 



(2.30) 



if the cut-off Rk is such that it diverges in this limit. If 
Rk does not diverge and is only very large. 



fc=A 



H . 



(2.31) 



2.2 The exact RG equation and 
its properties 

The RG equation on F^, sometimes called Wetterich's 
equation, that is the differential equation dk^k = fiXk) 
is derived in detail in the Appendix, section [3.3.21 The 
strategy is to obtain first an evolution equation for Zk^ 
then to deduce the equation on F^;. It writes 

dkVk = i / dkRkiq) (Fi'^[M] +7^fe) ' (2.32) 

where TZk{x,y) = Rk{x — y). The inverse 

(f^^^[A/] + TZk j has to be understood in the opera- 
^ ^ q-q 

tor sense. It is convenient in practice to rewrite Eq. (j2.32p 
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asb 



(2.33) 



where dk acts only on the /c-dependence of Rk and not 



on r 



(2) . 



dk 



dRk d 



dk dRk 

and the trace means integral over q (and for more com- 
plex theories, summation over any internal index). 



-q 



(2.34) Figure 2.3: Diagrammatic representation of the "full" 
propagator that is the M-dependent function Gk[M] = 



(ri^)[Af] 



2.2.1 Some general properties of the ef- 
fective average action method 

Let us mention some important properties of . 

• If the microscopic hamiltonian H and the functional 
measure are symmetric under a group G and if there 
exists a cut-off function Rk such that the mass term 
AHk respects this symmetry, then Tk is symmetric 
under G for all k and thus so is F = Ffe=o. It can 
happen that there is no mass-like term that respects 
the symmetry whereas the theory, that is F, is in- 
variant under G. This means that the symmetry 
is broken for all finite k and that the symmetry is 
recovered only for k ^ 0. This is the case of gauge 
symmetry. This symmetry breaking term can be 
controlled by modified Ward identities that become, 
in the limit k 0, the true Ward identities. It re- 
mains nevertheless difficult up to now to compute 
RG flows in a completely controlled way in gauge 
theories. 

• An exact RG equation for theories involving 
fermions can also be derived along the same line. 



Eq. ()2.33p looks very much like (the derivative of) a 
one-loop result since at one-loop: 

Ffc = i/+iTrlog(i/(2)+^^^ . (2.35) 

Thus, substituting i/^^) ^-j^g F^^'[Af] func- 
tional in the derivative with respect to k of Eq. (|2.35p 
changes the one-loop result into an exact one! There 
exists a diagrammatic representation of the RG 
equation written as in Eq. (|2.32[) and that empha- 
sizes its one- loop structure, see Figs l2.3l and l2.4l 




Figure 2.4: Diagrammatic representation of the RG 
equation on the effective average action. The line rep- 
resents the "full" propagator, Fig |2.31 the cross dkRk[q) 
and the loop, the integral over q. 



We define 

Gk[M] = (F^')[A/]+7^^ 



(2.36) 



which is the "full", that is Af-dependent, propagator, 
see Figl; 



• This one-loop structure has a very important prac- 
tical consequence: only one integral has to be com- 
puted and thanks to rotational invariance, it is one- 
dimensional. This is very different from pertur- 
bation theory where Z-loop diagrams require I d- 
dimensional integrals. This is a tremendous sim- 
plification of the NPRG method compared to per- 
turbation theory. 

• The perturbation expansion can be retrieved from 
the NPRG equation (and the all-order proof of 
renormalizability can be simpler in this formalism) . 

• Because of the dkRk term in the NPRG equation 
(|2.32p . only momenta q^ of the order of fc^ or less 
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contribute to the flow at scale k (we come back in de- 
tail to this point in the following). Thus the RG flow 
is regular both in the ultra-violet and in the infra- 
red. All the divergences of perturbation theory are 
avoided : we compute directly the RG flow and not 
first the relationship between bare and renormalized 
quantities from which is computed, in a second step, 
the RG flow. 

• k acts as an infrared regulator (for fc 0) somewhat 
similar to a box of finite size ~ /c~^. Thus, for fc > 0, 
there is no phase transition and thus no singularity 
in the free energy Tk- At finite k, everything is 
regular and can be power-expanded safely. We can 
therefore conclude that 

(i) the singularities of T build up as k is lowered and 
are thus smoothened by k in Ffe, 

(ii) the precursor of the critical behavior should al- 
ready show up at finite k for \q\ 3> k. 

• An important consequence of the regularity of 

at A: > is that it can be expanded in a power 
series of VM{x). For slowly varying fields M{x) 
this expansion is expected to be well-behaved. This 
is the basis of the derivative expansion that consists 
in proposing an ansatz for Tk involving only a finite 
number of derivatives of the field. [32] Two of the 
most used approximations, based on the derivative 
expansion, are 

Tk = j d'^x (^Uk{M{x)) + i (VM)'^ (2.37) 

called the local potential approximation (LPA) since 
no field renormalization in front of the derivative 
term is included and 

Tfe = y" d'^x (^Uk{M{x j) + ^Zk{M) (VM)'^ 

(2.38) 

often called the 0{d^) approximation or the leading 
approximation!! Of course, in principle, Tk involves 
all powers of VAf compatible with rotational and Z2 
invariance. We shall come back in detail on these 
approximations in the following. 

®The function Zj,(M) has, of course, nothing to do with the 
partition function Zf^[B\ introduced in Eg. 112.181 1 although it is 
customary to use the same name for both functions. 



• We shall see in the following that by working with 
dimensionless (and renormalized) quantities, the 
NPRG equation can be rewritten in a way that 
makes no explicit reference to the scales k and A. 
This will allow us to find fixed points from which 
quantities like critical exponents can be computed. 
The fact that A disappears from the flow equation 
also ensures that the "group law" of composition of 
RG transformations is satisfied. This self-similarity 
property of the flow will be automatic even if Ffc is 
truncated (as in Eqs. (|2.37l2.38p for instance). This 
is a major advantage of this method compared to 
many others where renormalizability, that is self- 
similarity, is spoiled by approximations (Schwinger- 
Dyson approximations for instance). 

• A property that follows from the last point is the 
decoupling of massive modes. Let us consider a 
theory with one very massive mode (mi) and an- 
other one with a lower mass (7712). Since only the 
modes |g| < k contribute to the flow of F^, once 
/c <C mi, the massive mode does (almost) no longer 
contribute to the flow : it decouples. This means 
that it contributes to the flow when k is between 
A and mi (or, in real space, when the running lat- 
tice spacing is between a and its correlation length 
rrii^) and almost not below toi. This also means 
that if we are given a model at a scale fco <C mi, we 
cannot know whether the underlying "fundamental 
theory" (at scales much larger than mi) involves or 
not a massive particle since below this scale there 
exists almost no signal of the presence of this mode 
in the theory (that is, in the RG flow). Except for 
very precise measurements, the high energy sector of 
the theory has been decoupled from the low energy 
one and its contribution to the theory at low energy 
amounts to renormalization of the couplings occur- 
ing during the part of the RG flow corresponding 
roughly to mi < fc < A. 

• Self-similarity and decoupling of massive modes will 
have universality (in the vicinity of a second order 
transition) as a consequence. We shall see that uni- 
versality is not just the consequence of the existence 
of a fixed point. It is the consequence of a very pecu- 
liar geometry of the RG flow, see section \2l6\ [49l [5T] 
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2.3 Approximation procedures 



The NPRG equation (|2.32p is an extremely complicated 
equation. It is a functional partial differential equa- 
tion since it involves the functionals Ffc [Af ] and F^^'' [M] . 
Needless to say that we do not know how to solve it in 
general. Some approximations are thus required. Two 
main kinds of approximations are usually considered : 
the Green function approach and the derivative expan- 
sion. In both cases, the strategy consists in solving the 
RG equation in a restricted functional space and not as 
a series expansion in a small parameter. This is why we 
can hope to obtain non-perturbative results. Both meth- 
ods need to project consistently the exact RG equation in 
the functional space that has been chosen. Of course, the 
quality of the result will depend crucially on the choice 
of space in which we search for a solution. Depending on 
the problem, one choice can be drastically better than 
another. In all cases, it is impossible to know whether we 
have missed some physically crucial ingredient by mak- 
ing one choice rather than another one. But this problem 
is not specific to the NPRG. It is generic in physics. . . 

Let us review briefly the Green function approach and 
then explain in some details the derivative expansion 
which is the most employed method in statistical me- 
chanics. 



Therefore, setting Rk = dtRk = kdkRk we obtain 



5Mp 6Mp> 



Rk,q^ 



.(3) 



^k q^.-q^'^ k p' ^q^^-q^'^k qr^^-q-i 



9 / ^k^qi Gk qi,^q2^k q^-qs.p.p'^^^ls-qi 
^ ■■'{9.} 



(2.41) 



where both the left and the right hand sides are func- 
tions of Mq since they have not yet been evaluated in a 
particular configuration. These equations look terrible 
but in fact they are not since there exists a diagram- 
matic way to obtain them automatically. We represent 
rl^^ as in FigHH 

k pi,p2,P3 ° 



n(3) 




P3 



Figure 2.5: 

r?^ [Ml. 



Diagrammatic representation of 



We obtain for StF^^-* the graph in Fig, 



2.3.1 The Green function approach 



From Eq. (|2.32p we can deduce the infinite hierarchy of 
RG equations on the correlation functions defined by 



=,(n) 
k,pi...p„ 



<5"Ffc 



SM„ 



. SAL 



(2.39) 




taken in a particular field configuration (the zero and the 
uniform field configurations being the most useful ones). 
We define t — log A;/ A which is often called the RG 
"time". From Eq. (l2.32|) . we obtain 



6Tk 

SMr, 



Rk{qi)Gk 



5T 



(2) 



92:-'?3 



«l:-92 



Gk , 



(2.40) 



Figure 2.6: Diagrammatic representation of (?tF[,^^. 

Of course, if we evaluate dtT^j!^^ in a uniform field con- 
figuration, the momentum is conserved at each vertex 
and for each "propagator". The function is thus non- 
vanishing only at zero momentum and we obtain the 
equation in Fig l2.7l 

It is clear on the diagrammatic representation that 
only one momentum integral remains. The RG equation 
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Figure 2.7: Diagrammatic representation of dtTj^ eval- 
uated in a uniform field configuration. 



for dtTj^ evaluated in a uniform field configuration is 
given in Fig |2.8l 




Figure 2.8: Diagrammatic representation of dtT^. eval- 
uated in a uniform field configuration. Note that the 
factor i has not been represented on the figure for the 
tadpole because it can be retrieved from the topology 
of the graph itself (in fact the minus sign can also be 
retrieved from the graph). 

It is clear from these diagrammatic representations 
that dtV'-^^ involves r["+^^ and r^"+^\ If we want to 
solve this infinite tower of equations, we have to trun- 
cate it. A possible truncation consists, for instance, in 
keeping r^^-* and r^"*-* evaluated at Af = and to ne- 
glect the contribution of F^^-* in the equation on ^tF^,"*-* 
(rf^[M = 0] = in a Z2-invariant model). A better 
method is to find an ansatz for F^^"* in terms of F^^-* and 
F^"*-*. In both cases, the systems of equations become 
closed and can, at least in principle, be solved. Let us 
finally notice that this method consists in truncating the 
field-dependence of Tk [M] while keeping the momentum 
dependence of F^^'' and F^''-' . Tremendous improvements 
of this type of approximation has been performed these 
last years. [51 [51 [Ml 



We now study another truncation method which is 
somewhat the reverse. 

2.3.2 The derivative expansion 

The principle of this approximation has already been 
introduced previously, see Eas. (|2.37l2.38p . The underly- 
ing idea is that we are mostly interested (for the study 
of critical phenomena) in the long distance physics, that 
is the |g| —i- region of the correlation functionsl^ Thus, 
we keep only the lowest orders of the expansion of F^ in 
VM while we keep all orders in the field M 

Ffe = y" d'^x (uk{M{x)) + ^Zk{M) (VAf )'^ + 0(V4) 

(2.42) 

This approximation is based on a somewhat opposite 
philosophy as the one that prevailed in the Green func- 
tion approach. However, it should be clear that for sta- 
tistical mechanics, the most important informations — 
e.g. the equation of state — are hidden in the effective 
potential Uk=o of the theory that, therefore, needs to be 
computed as accurately as possible (see however [53|). 

It is in fact remarkable that we can combine both 
methods by making a field expansion of C/fc, Z^, . . . on 
top of the derivative expansion while preserving many 
non-trivial results. The simplest such truncation con- 
sists in using the LPA and in keeping only the first two 
terms of the expanion of Uk in powers of 

Tk[M] = J d''x{g2,k + 54,fc + ^(VM)^) (2.43) 

With this kind of ansatz, the RG equation on F^ becomes 
a set a ordinary differential equations for the couplings 
retained in the ansatz: 

dtdn.k = f3n i{gp,k}) ■ (2.44) 

^The computation of quantities like the total magnetization 
or the susceptibilities require only the knowledge of the spin-spin 
correlation function at zero momentum. The same thing holds for 
the correlation length. 

^''Let us notice that if the fc-dependence of the couplings was ne- 
glected, this ansatz would exactly coincide with the ansatz chosen 
by Landau to study second order phase transitions. We know that 
it would lead to the mean field approximation. It is remarkable 
that keeping the scale dependence of the couplings and substitut- 
ing precisely this ansatz into the RG equation of Fj, is sufficient to 
capture almost all the qualitative features of the critical physics of 
the Ising and 0(N) models in all dimensions (see the following). 
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For instance, if the truncation above is considered we 
fincS 

dtKk = -id^2)Kk + 6vdlf{2Kk\k) (2.45) 

(2.46) 

dtXk = (d- 4)Afe + 6vd Xl lli2KkXk) (2.47) 

where if, I2 are defined in the Appendix, section 13.11 
These equations are aheady non perturbative since the 
functions , I2 are non polynomial. 

If t/fe, . . . are not truncated in a field expansion, 
the RG equation on becomes a set of coupled par- 
tial differential equations for these functions. The ini- 
tial condition at scale A is given by the hamiltonian of 
the model. Before studying in detail this approximation 
method let us make an important remark. 

Contrarily to perturbation theory where only the 
renormalizable couplings are retained in the renormal- 
ized action, all powers of the fields appear in the ansatz 
of Eq. (|2.42p . There is no longer any distinction — at 
this level at least — between the two kinds of couplings, 
renormalizable and non renormalizable. This point to- 
gether with a discussion of the notion of perturbative 
and nonperturbative renormalizability will be discussed 
in the following. 



Eq.(I131]). The only "difficulty" of this calculation is to 
invert F^ [M]+Rk for uniform field configurations. This 
is where truncations of F^ are useful: they enable explicit 
calculations. For the FPA ansatz this is performed in 
detail in the Appendix, section 13.3.31 The final result 
reads: 



dkUkip) = 



where 



1 



2 g2 



dkRkiq) 



Rk{q) + Ul{p) + 2pU'ap) 



1 



(2.49) 



(2.50) 



is the Z2-invariant and Ul.{p) and U\^{p) are derivatives 
of Uk with respect to p. An important remark is in order 
here. Once the angular integral has been performed, the 
integrand of Eq. ()2.49|) is of the form (up to a constant 
factor) 



f{q\w) 



dkRkiq^ 



w 



q- + Rkiq^) 
Fet us consider a typical cut-off function 

„2 



Rkiq' 



1 



(2.51) 



(2.52) 



Then, for generic values of w, the typical shape of / at 
fixed k is given in Fig l2.9l (for d > I). 



2.4 The local potential approxi- 
mation for the Ising model 

2.4.1 The flow equation of the potential 



We now consider the ansateEa. (|2.37p for a Z2-symmetric 
theory. As already mentioned, the problem is to project 
the RG equation (|2.32p on the potential Uk- This is natu- 
rally performed by defining the potential as F^ computed 
for uniform field configurations: 



C/fe(M, 



unif. , 



n 



Ffe[Mu„if.] 



(2.48) 



where is the volume of the system. To compute the 
RG flow of Uk we act on both sides of this equation 
with dt and we evaluate the right hand side thanks to 




^'^The derivation of these equations is outlined in section [2IZ21 



Figure 2.9: The typical shape of the integrand f{q)- 

As expected, only a window of momenta around k con- 
tributes to the flow at scale k. We see in particular that 
the rapid modes are efficiently integrated out by this kind 
of cut-off function. This explains the decoupling of mas- 
sive modes already discussed section 12.2.11 The cut-off 
function (|2.52p has been used many times in the litera- 
ture but it turns out that another one is sometimes more 
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convenient because it allows to perform analytically the We define the dimcnsionless variables by 
integral in Eq. (|2.49[) . It writes 



Rk{q^)^{P -q^)e{P -q^) 



(2.53) 



where 9 is the step function. With this choice of Rk we 
easily find 



Ul{p) + 2pU'^{p) 



(2.54) 



where Vd is a numerical factor defined in Appendix, sec- 
tion [SH] All we can learn about the model at this ap- 
proximation is contained in the solution of this equation. 
We can already see that, as expected, this equation does 
not admit a fixed point potential U* . We have already 
seen that it is necessary to go first to the dimcnsionless 
variables to find a fixed point. This is what we now do. 



2.4.2 The scaling form of the RG equa- 
tion of the dimensionless potential 

We have already emphasized when studying block-spins 
that it is necessary to measure all lengths in units of the 
running lattice spacing to find a fixed point of the RG 
flow. Looking for a fixed point is a convenient method 
to study the critical behavior of a model and we shall 
now spend some time deriving the RG equation on the 
dimensionless potential. 

In our formalism, k is the analogue of the inverse 
running lattice spacing and, to find a fixed point, we 
must therefore "de-dimension" all dimensionful quanti- 
ties thanks to k. This is equivalent to measuring all 
lengths in units of the running "lattice spacing" k~^ . We 
have 



[Tk] = fc° 



[M] = fc"^ and [Uk\ = k'^ 



q^ 




(2.56) 




(y) = k^yrijj) 


(2.57) 


X = kx 




(2.58) 


M{x) = k^ 


M{x) 


(2.59) 


iJt{M{x)) - 


k-''Uk{M{x)) . 


(2.60) 



To derive the RG equation on Ut we must keep in mind 
that the derivative dt in Eq. (|2.54p is computed at fixed p 
whereas we want now to compute it at fixed p. The de- 
tailed calculation is performed in the Appendix, section 
13.3.31 and the result is: 



dtUt ^-dUt + {d-2)plJ[ 
-2vd / dyy'^/^+i 



r'{y) 



y{l + r{y)) 



2pUi' 



(.2.61) 



Once again, with the cut-off (|2.53p the integral can be 
performed analytically and we find 



dtUt - -dUt + {d^2)pU^ + 



4wd 



1 



d l + U[ + 2pU[' 



(2.62) 



(2.55) 



We clearly see on this equation that the fiow of Ut has 
two parts, one that comes from the dimensions of Uk and 
p and one that comes from the dynamics of the model. 
This RG equation on Ut is a rather simple partial differ- 
ential equation that can be easily integrated numerically. 
We are thus in a position to discuss the critical behavior 
of the Ising model and to look for fixed points. [32] 

2.5 The critical and non-critical 
behavior of the Ising model 
within the LPA 

Having derived the RG flow of the potential, we can 
relate a "microscopic" model defined at scale A by an 
hamiltonian (or directly by Fa) with the free energy T = 
Tfc^o • Let us emphasize that there is no reason why the 
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hamiltonian H should involve only and (j)'^ terms and 
not (ffi, (j>^, . . . terms. In fact, the Hubbard- Stratonovich 
transformation enables to obtain the exact potential at 
the scale of the lattice spacing of a magnetic system. In 
the Ising case, the potential thus obtained is 

Ua{(I>) oc log cosh (/) (2.63) 

and is thus non-polynomial. In the NPRG framework 
this is not a problem since, even if J7a(0) was a polyno- 
mial at scale A, it would become non-polynomial at any 
other scale (all Z2-invariant couplings are generated by 
RG transformations, see section \2M for a discussion). 

Let us anyway, for the sake of simplicity, consider a 
dimensionless potential at scale A of the formP^ 

UA{p)^^{p-^Af (2.64) 

with KA > 0. At the mean-field level and for uniform 
fields, = and we would deduce at this approx- 
imation from Eq. (|2.64p that the system is in its broken 
phase with a spontaneous magnetization per unit vol- 
ume Msp — A^a/2kaF^ However, the integration of 
the fluctuations can drastically change this picture: the 
minimum of the potential has a non-trivialRG flow that 
can drive it to zero. If this occurs, the system is in fact 
in its symmetric (high temperature) phase. 

Let us call K{k) the running minimum of the di- 
mensionless potential at scale k (more appropriately at 
"time" t = logfc/A): 

5pt/feu, =0 (2.65) 

It is physically clear, and this can be checked on the 
flow of f/fe, that the spontaneous magnetization decreases 
in comparison with its mean-field value because of the 
fluctuations. This means that the true spontaneous mag- 
netization is always less than the mean-field spontaneous 
magnetization ^/2ka. There are thus three possi- 
bilities. 

^^The two models corresponding to the two initial conditions 
Eq |2.63l and Ea l2.64l belong to the same universality class. Thus 
they both have the same set of critical exponents. However, they 
differ as for their non universal quantities. 

^^At vanishing external magnetic field, the magnetization is 
given by -B = = dU / dM and corresponds therefore to the lo- 
cation of the minimum of the potential. 



(i) The system is in its broken phase (low temperature) 
and the (dimensionful) spontaneous magnetization den- 
sity is given by 

M,p = V^poik = 0) (2.66) 

where po{k) is the location of the minimum of the (di- 
mensionful) potential C/fe(/9)Ill The relation between 
K,{k) and po{k) is 

po(fc) = k'^-^K{k) (2.67) 

from Eq. (12311). Thus, if po{k) Msp/2 when fc -> 0, 
K(fc) has to flow to infinity (for d > 2) as k"^^'^ in this 
limit. 

(ii) The system is in the high temperature phase and 
the spontaneous magnetization is vanishing. Thus, as 
k decreases the minimum K{k) must decrease and, at a 
flnite scale fcg, must hit the origin. It is easy to guess 
that ks must be of the order of the inverse correlation 
length since as long as <C ^ the coarse-grained sys- 
tem remains strongly correlated and still behaves as if it 
was critical. It is only when k~^ ^ ^ that "the system 
can feel" that its correlation length is finite and that its 
magnetization is vanishing. 

(iii) The system is critical. The spontaneous magnetiza- 
tion is vanishing which means that po{k) — k'^~'^K{k) — > 
as fc 0. Note that this does not require that k(0) = 
since K,{k) is multiplied by a positive power of k, at least 
for d > 2. In fact, K.{k) reaches a finite fixed point value 

K*. 

These three cases are summarized in Fig l2.10l 
For a generic initial potential 

C7a(M) = ^{~p - + ^{p - n^f -t- . . . (2.68) 

the critical value of ka is a function of all the other 
couplings Aa,M3,a, ■ • ■ This means that in the space of 
dimensionless coupling constants there exists a "critical 
hypersurface" of co-dimension 1 that corresponds to sys- 
tems that are critical. Generically, ka is a regular func- 
tion of the temperature around Tc and at first order we 
can assume that 

KA ~ Tc — T . (2.69) 

-'^''We shall see in the following that there is a subtlety here be- 
cause of the convexity of the potential in the limit fc — > 0. 
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Figure 2.10: Behavior of the running minimum ^(fc) 
of the dimensionless potential. From top to bottom: 
T < Tc, T — Tc and T > Tc- Kc is the critical initial 
value of K\ for which the system is critical. This value 
should not be confused with k* which is the fixed point 
value of K{k). In the low temperature phase the dimen- 
sionless running minimum diverges whereas the dimen- 
sionful minimum converges to the value of the sponta- 
neous magnetization since it is multiplied by a positive 
power of the scale k that compensates exactly the diver- 
gence of K{k). 

This allows us to relate the coefficients of the microscopic 
hamiltonian to the temperature. Note that contrarily to 
the mean-field analysis for which criticality is reached 
when the coefficient tq of the quadratic term of the po- 
tential (that is ro0^/2) is vanishing: ro cx T — Tc, this is 
not true here since criticality does not correspond to the 
vanishing of the bare mass term. The mean-field analysis 
is wrong in this respect. 

Let us now show what we expect for U = Uk=o and 
for iJk=o- 

2.5.1 The low and the high temperature 
phases 

In the low temperature phase, we expect a spontaneous 
magnetization, either up or down. More precisely, the 
equation 

is expected to have a solution +Msp for i? ^ 0+ and 
—Mgp for B 0^ . Moreover, U must be a convex func- 
tion of M since it is obtained from a Legendre transform 



of W which is convex. Thus, U must have a very peculiar 
shape since it must have two minima at ±Msp and no 
maximum in between (otherwise it should not be con- 
vex). The only possibility is that it is flat in between 
-Msp and +Msp. 



u 




Figure 2.11: Shape of the effective potential U = Uk=o 
for T < Tc. 

The convexity of the effective potential is not pre- 
served by perturbation theory: the convex envelop has 
to be taken by hand. In fact, it is notoriously difficult to 
compute "safely" a convex effective potential and this is 
a property that is reproduced by the NPRG already at 
the LPA. Note that this is no longer the case if a field ex- 
pansion of the potential is performed (a polynomial can 
never be flat on a whole interval). In fact, starting with 
a potential Ua showing a double- well structure and with 
parameters ka, Aa, . . . such that the system is in its low 
temperature phase, the RG flow of Uk is such that as k 
is lowered [32] 

• the minimum po{k) moves towards the origin and 
eventually reaches a limit equal to the spontaneous 
magnetization 

• the maximum of the potential located between the 
two minima decreases and goes to as fc (the 
"inner part" of the potential flattens). At fc = 0, 
Uk=o looks like the potential of Fig l2.11l 

As T —!■ Tc, that is as k\ — > Kc, Msp moves towards 
the origin and aX T — Tc coincides exactly (by deflniton 
of Tc) with the origin. In the high temperature phase, 
that is KA < Kc, the spontaneous magnetization van- 
ishes. Thus U has a single minimum at the origin in this 
case. 
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Let us finally notice that when k ^ 0, Uk is not convex 
in general. This is normal since Tk is not the Legendre 
transform of Wk {Wk is convex) because of the term 
1/2 J RkMM that has been subtracted, Eq.^^. 

2.5.2 The critical point 

For ka — the minimum of the potential Uk (or Uk) 
for fc > is non vanishing. It is only for fc — > that the 
minimum po{k) of Uk reaches the origin. This means 
that the minimum of the potential never stops running 
ai T — Tc and that fluctuations of all wavelengths are 
required to make po(fc) vanish. 

To characterize the critical point it is interesting to 
work with Uk instead of Uk since, at this point, the long 
distance physics (compared with A~^) is scaleless. This 
means that the potential, properly rescaled thanks to fc, 
should be fc-independent at T = a-nd for sufflciently 
small fc: it must be a fixed potential U*{p) 



U 



dtU*{p) = 0. 



(2.71) 



We deduce from this equation and from Eqs. (|2.53|2.62p 
that U* (p) is solution of 



= -dU* + {d-2)pU* 



1 



d 1 + U*' + 2pU*" 

(2.72) 

At first sight, the situation looks paradoxical since this 
is a second order differential equation that should ad- 
mit infinitely many solutions indexed by two numbers 
whereas we expect only one fixed point in d = 3 cor- 
responding to the universality class of the Ising model. 
In fact, it can be shown that among all these solutions, 
only one is well defined for all p S [0, oo[.[T3 

An easy way to show this is the "shooting method". 
We first write down the fixed point equation for the 
derivative of U* with respect to M: U* '{M) analogous 
to Eq. (|2.72p and which is more convenient. Then, one 
inital condition is given by the Z2 symmetry: U* '{M = 
0) = 0. The other initial condition: U*"{M = 0) is 
then adjusted so that U* '(M) exists for all M. For a 
generic U* "(M — 0) this is not the case: at finite M, 
U* ' either blows up at -t-oo or at —00. By dichotomy, 
the value U* "{M = 0) can be fine-tuned so that this oc- 
curs for larger and larger M . The fixed point potential 
thus obtained is shown in Fig l2.12l 




Figure 2.12: Dimensionless fixed point potential of the 
Ising model in d = 3. The minima of this potential are 
located at ±'\/2k*. The spontaneous magnetization is 
vanishing although 7^ (see the text). 



The same fixed point potential can be found "dynami- 
cally" by integrating the RG flow and by fine-tuning the 
value of KA to get closer and closer to Kc- We find that for 
an initial ka very close to Kc, the potential Ut spends a 
very long RG time close to U* before either departing in 
the high or low temperature phase. Thus, all dimension- 
less couplings reach a plateau before blowing up. On this 
plateau, we obtain a very good approximation of U*{M) 
very close to the one found by the shooting method. For 
K{k) very close to Kc this is represented in Fig |2.13l 



Figure 2.13: Behavior of K(fc) very close to Kc- 

One can observe that during a transient regime, K{k) 
for fc ~ A is not stationary even at the critical tempera- 
ture. This regime simply corresponds to the RG "time" 
necessary to approach the fixed point. It is non- universal 
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since it depends on the starting point Kc on the critical 
surface. 

Let us finally point out a subtlety. When A; > 0, it is 
possible to reconstruct Uk{M) from Uk{M) by a some- 
what trivial rescaling: 

Uk{M) = k''Uk(k-^M^ . (2.73) 

The k factor acts as a magnification scale when we go 
from M to M (at least for d > 2). For fc ^ an in- 
finitesimal range of M around the origin is mapped onto 
a finite range of M. It is therefore possible — and this is 
what indeed occurs at T = Tc — that Uk=o has a trivial 
shape around M = showing a minimum only at (see 
Fig l2.11l for Msp — > 0) whereas Uk=o shows a double-well 
structure (see Fig l2.12p ! 



2.5.3 The critical exponents 

Within the LPA, there is no possibility of modifying the 
mean-field like g^-dependence of T''^\q) (evaluated at 
M = 0) since there is no renormalization of the deriva- 
tive term. The exponent rj is therefore vanishing at 
this order of the derivative expansion (see Eq. (|1.70p and 

r(2)(g) = i/G(2)(g)). 

r/LP^ = . (2.74) 



The exponent /3 defined in Eq. (|1.85p can be calculated 
directly by fitting the behavior of Mgp defined in Fig l2.11l 
as a function of ka — which is itself proportional to 
Tc — T. Of course, this is not very convenient although 
feasible. 

As for the exponent we have seen that it is related 
to the behavior of the RG flow around the fixed point, 
Ea. (|1.74p and, more precisely, is the inverse of the posi- 
tive eigenvalue of the flow at the fixed point. This means 
that very close to the fixed point and away from the crit- 



ical surface, the potential is such thalF^ 

U{M, t) ^ U* (M) + e e^*/''it(M) . (2.75) 

It can be shown that u{M) must behave as a power law 
at large M and that this ensures that there is a unique 
value of ly such that (P?75)l holds. With the ^-cut-off 
function, Ea. (P35)) . one finds at d = 3: ^21^1 ^5] 

i/LP^ = 0.65 (2.76) 

to be compared with — 0.6297(5) obtained by Monte 
Carlo simulations. 

Several remarks are in order here. 

• In principle, from v and rj all the other critical ex- 
ponents can be calculated thanks to the scaling re- 
lations derived in section [1.5.31 It is interesting to 
verify that these relations are not spoiled by the 
LPA. Indeed, by computing separately all the expo- 
nents, it is found that the scaling relations among 
them are very well satisfied. 

• In the exact theory, no physical result should de- 
pend on the choice of cut-off function Rk- However, 
once approximations are performed a spurious de- 
pendence on the choice of Rk is observed. As a 
consequence, the whole scheme makes sense only if 
this dependence is weak. It is of course very difh- 
cult to obtain general results on this point since the 
space of cut-off functions is infinite dimensional and 
that we cannot sample it efficiently. In practice, we 
should choose a space of "reasonable" cut-off func- 
tions and study the variations of physical quantities 
like critical exponents in this space. This has been 
done in some details and it is observed that the de- 
pendence of ly on Rk is indeed weak. [331 351 135] Let 
us emphasize that this problem is not specific to 
the NPRG. Even in the perturbative schemes the 
critical exponents acquire a spurious dependence on 
several choices made during their calculations. 

^^The following relation is not general since if we choose a point 
on the critical surface, the corresponding potential is attracted to- 
wards U*{M). This approach is governed by the so-called critical 
exponent w corresponding to correction to scaling. Thus, in gen- 
eral, infinitesimally close to the fixed point, the evolution of the 
potential is given by the sum of two terms, one describing the ap- 
proach to U* on the critical surface and one describing the way 
the RG flow escapes the fixed point if one starts close but away 
from the critical surface. 
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• The LPA is certainly not appropriate for the com- 
putation of critical exponents in c? = 2. We have 
seen several times that this dimension plays a spe- 
cial role in the formalism although nothing spectac- 
ular is expected in d = 2 for the critical behavior, 
at least for the Ising model. This comes from the 
fact that 77 = starts to be a bad approximation 
at and below d = 2 for this model (it is worse for 
the 0{N) models with N > 2 because of Mermin- 
Wagner's theorem). The exact value of 77 is known 
in c? = 2 from Onsager's solution of the Ising model: 
r] = 0.25. At low dimensions, going to the next order 
of the derivative expansion cures most of the prob- 
lems encountered at the level of the LPA. [32l [56] 

• As explained in the following for the 0{N) models, 
it is possible to go beyond the LPA and to compute 
with greater accuracy the critical exponents, 77 in 
particular. This kind of calculations has been per- 
formed by several authors [57] at order 0{d^) of the 
derivative expansion for the Ising and 0(N) models 
in d = 3 and d = 2 [58] and also at 0(9^) for the 
Ising model in d = 3. Let us quote the best results 
obtained for the Ising model in d = 3: 



order 




V 




0.6506 







0.6281 


0.044 




0.632 


0.033 


7-loops 


0.6304(13) 


0.0335(25) 



Table 2.1: Critical exponents of the three dimensional 
Ising model. 9", and 9^ correspond to the order of 
the truncation of the derivative expansion (the NPRG 
method) [46]. For completeness, we have recalled in the 
last line, the results obtained perturbatively[5]. 



Let us now review some conceptual aspects of renor- 
malization and in particular what Wilson's RG teaches 
us about peturbative and non perturabtive renormaliz- 
ability. 



2.6 Perturbative renormalizabil- 
ity, RG flows, continuum 
limit, asymptotic freedom 
and all that. . . 

Our understanding of the physical meaning of renormal- 
izability has drastically changed with the works of Wil- 
son and Polchinski[T2] (and many others)0. It is even 
astonishing to see that while the technical aspects of 
perturbative renormalization have not changed since the 
50 's, the interpretation of the notion of renormalizability 
has been deeply modified these last twenty years. Let us 
mention just a few points that illustrate this change of 
viewpoint on renormalization: 

(i) The regularization scale of a field theory — the cut- 
off of the Feynman integrals or the inverse lattice spacing 
for instance — was considered as an unphysical scale in- 
troduced during intermediate calculations and that had 
to be sent to infinity at the end of the calculation. It 
is now thought to be the scale where new physics takes 
place and is thus a fundamental scale of the model. Most 
of the time, it is not necessary to send this scale to in- 
finity (but perhaps for technical convenience). 

(ii) For the same reason the "bare theory" (associated 
with the scale of the cut-off) was considered as unphysi- 
cal whereas the renormalized theory was supposed to be 
fundamental. The bare theory is now often considered as 
the fundamental theory whereas the renormalized theory 
is only an effective theory valid at long-distance (long- 
distance meaning distances that are large compared to 
the inverse ultra-violet (UV) cut-off). 

(iii) The only valid theories were thought to be built 
with renormalizable couplings only while nonrenormal- 
izable couplings were believed to be forbidden (for rea- 
sons that looked rather mysterious). Nonrenormalizable 
couplings were thus considered as terrible monsters that 
nobody was able to master while super-renormalizable 
couplings were considered as very kind beasts. In this 
zoo, renormalizable couplings were classified just in be- 
tween these two kinds of animals. We now know that it 
is useless to incorporate nonrenormalizable couplings in 

^^Much of what follows in this section comes from the works of 
Bagnuls and Bervillier to whom I owe much of my understanding 
of this subject. 1151 [ST) 
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the "bare" action — that is in the UV regime — at least 
if one is interested only in the "universal" behavior of the 
model under study (see the following for a precise expla- 
nation of this statement). Put it differently, the contri- 
butions of the nonrenormalizable couplings included in 
the bare action die out at large distances and are thus 
under control. Reciprocally, super-renormalizable cou- 
plings of the bare action (when they exist) do not make 
any problem in the ultra-violet regime but are those that 
dominate the large distance physics of the theory and are 
thus the most important terms to control in this regime. 

(iv) Removing the cut-off, that is taking the "infinite 
cut-off limit" order by order of perturbation expansions, 
was supposed to be the real issue of renormalizability. 
It turns out that the possibility of taking the infinite 
cut-off limit at all orders is not equivalent in general to 
the possibility of taking the "continuum limit", that is of 
defining the field theory at the limit where the nonper- 
turbative regulator (the lattice for instance) is removed. 
Moreover, taking the continuum limit of a model is not 
in general a physically relevant issue since at very short 
distances the relevant physical model can well be some- 
thing else than a field theory. This is already the case 
for the Ising model. The fundamental theory can also be 
a field theory but involving many other terms than the 
renormalizable ones and even infinitely many. The cj)'^ 
model for instance is a good effective theory of the Ising 
model at large distances but the continuum limit of this 
field theory teaches us nothing about the small distance 
behavior of the Ising model. 



Let us now discuss in more details these structural 
aspects of renormalization. As a starting point let us 
consider the different status of the couplings in the per- 
turbative and nonperturbative frameworks. In the for- 
mer some couplings are called renormalizable while the 
others are called nonrenormalizable whereas in the latter 
they are treated on the same footing. We shall mainly 
study the massless case in what follows although it is 
possible to generalize almost all arguments to the mas- 
sive case. 



2.6.1 Renormalizable couplings and 
large river effect 

To understand the origin of the difference between renor- 
malizable and nonrenormalizable couplings, it is neces- 
sary to remember that in perturbative field theory, the 
action (bare or renormalized) is not supposed to be a 
physical quantity. It is only the "mathematical tool" that 
generates the vertices and the propagators of the theory 
from which Feynman diagrams are computed. As for the 
Feynman diagrams they are the building blocks of the 
Green functions that are physical (in particle physics, 
the 5'-matrix is physical). The effective action r[M] is 
also physical. The action S or the hamiltonian H of per- 
turbative field theory can well be polynomial and involve 
only a few couplings, T[AI] always involve all powers of 
the field in a non-trivial way. In the particle physics 
language and for the 0^ theory, this means that the con- 
nected Green functions G^"-* with n — 6,8,..., corre- 
sponding to the scattering of n particles are non vanish- 
ing and cannot be factorized into products of G^^-* with 
p < n corresponding to the scattering of fewer particles. 
This is a signature of the fact that a field theory involves 
infinitely many degrees of freedom that require infinitely 
many independent correlation functions to be faithfully 
describedl^ Thus all couplings gn — defined as the val- 
ues of the r("''s for some configuration of the external 
momenta — are non-trivial, even those corresponding 
to an arbitrarily large n, that is to an arbitrarily large 



^'^At first sight, this statement could seem incorrect in particle 
physics if one considers a given diffusion process. For instance 
the reaction e+e~ — > 77 seems to involve four bodies only. This 
is actually wrong since, as virtual states involved in the loop ex- 
pansion, an arbitrarily large number of particles can be exchanged 
during this reaction. The full Fock space structure is thus nec- 
essary to describe any kind of diffusion in the quantum and rela- 
tivistic framework. The same is true in statistical mechanics. The 
infinite number of degrees of freedom of the system is not the rel- 
evant point. A perfect gas for instance can well involve infinitely 
many degrees of freedom, we all know that the whole machinery 
of field theory is not necessary to study it. The important point is 
the number of degrees of freedom that effectively interact together, 
that is the value of the correlation length. As long as ^ ~ a field 
theory is not necessary since the system breaks down in small sub- 
systems of size ^ that are almost independent of each other. This 
is the reason why the law of large numbers is valid in this case 
and the fiuctuations gaussian. Field theory is relevant only when 
g 3> a in which case, for length scales I such that a <C i <5C ^ field 
theory is relevant. 
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Figure 2.14: Flow on the critical surface. The line i, 
called the large river, is an attractive submanifold of the 
RG flow of dimension one (in the Ising case) . It starts at 
the Gaussian fixed point and ends at the Wilson-Fisher 
fixed point (WF). All RG trajectories approach L very 
rapidly while the flow on L is "slow". All the starting 
points of the RG trajectories correspond to different ini- 
tial conditions, that is, at microscopic scales, to different 
Z2-invariant systems. The part of the RG trajectories 
that is far away from L is non-universal and corresponds 
to a transient regime. The axis {17^} represents the infi- 
nite number of axes different from 34 which is the cou- 
pling of the (j)^ term. The axis m represents the direction 
of the mass which is a relevant direction of the RG flow. 

number of fields. 

In Wilson's framework, all couplings gn have a non- 
trivial RG evolution that, a priori, needs to be taken into 
account. This is the reason why, for instance, we must 
keep a complete function of M for the potential Uk and 
not just the first terms of its expansion in M, that is the 
A'P and terms^l If at scale A we had retained only 
these couplings in the initial condition of the RG flow on 
Uk, we would have anyway found that all others would 
have been generated at any scale k lower than Thus, 

^® Remember that Uk corresponds to the zero momentum con- 
figuration of all r^"' . 

Truncating the field dependence of the running potential Uk 
by keeping only the Af^ and terms at all scales k is nothing but 
a very crude approximation that eventually leads to neglecting all 
functions r(") with n > 6 in the efl'ective action r[A/] = Fk^olM]. 



at first sight, all couplings should be treated on the same 
footing and the distinction between renormalizable and 
nonrenormalizable couplings seems to have disappeared 
in Wilson's approach. 

This is paradoxical because this seems to be in con- 
flict with what we know from perturbation theory. In 
the perturbative scheme, all F'^"^ are also nontrivial and 
thus all (renormalized) couplings corresponding to the 
values of the renormalized F^^ (at zero momentum for 
instance) are nontrivial. However, since all renormalized 

functions F^'^ can be expanded in powers of the renor- 
malized coupling (?4,/j, it is possible to compute all renor- 
malized couplings g6,R, ga,R, etc in terms of this unique 
coupling constant. This is what makes a renormalizable 
field theory predictive: once a finite number of (renor- 
malized) masses and couplings have been determined by 
external means (e.g. experiments) the whole theory is 
entirely determined and infinitely many quantities can 
be computed out of it. 

The paradox is resolved by explicitly computing the 
RG flows of the couplings gn,foT instance within the 
local potential approximation^ At sufficiently long dis- 
tance, that is for k sufficiently small, all RG trajectories 
— that take place within the infinite dimensional space 
of coupling constants — are found to be attracted to- 
wards a submanifold of dimension two, see Fie: l2.14l for 
the critical case where the attractor is of dimension one 
since the massis eliminated!^ 

The RG time necessary to reach this submanifold is 
very short whereas the flow within it is slow. This is 
particularly clear, and has been studied in detail by Bag- 
nuls and Bervillier, for the Z2-invariant theories in rf = 3 
when they are critical. [49l [5T] In this case, the RG tra- 
jectories belong to the critical surface. All of them, after 
a transient regime, almost collapse on a line, that we 
call L, joining the gaussian fixed point to the non-trivial 
fixed point describing the phase transition of the Ising 
model (called the Wilson-Fisher fixed point). Thus, in 
the long-distance physics compared to A^^, that is be- 

All qualitative features that are explained below do not depend 
on this approximation. 

^^In general, the attractive submanifold is of dimension the 
number of (perturbatively) renormalizable couplings including the 
masses. In the critical case, the mass is vanishing by definition 
and the attractor is of dimension one for models for which only 
one coupling is renormalizable. 
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yond the transient regime, the RG flow behaves as if it 
was driven by a unique "couphng". As long as L has a 
non singular projection on the axis (74 corresponding to 
the M''-coupling, it is possible to describe the flow along 
the line L from the flow of the (74-coupling alonel3 This 
is what perturbation theory does. This explains why in 
this framework all couplings have on one hand a non- 
trivial RG flow and on the other hand a flow determined 
by (?4 alone. 

Bagnuls and Bervillier have used the following 
metaphor: [49l [51] the RG trajectories on the critical sur- 
face are like rivers in the mountains. In the valley, there 
is a large river L along which the flow is slow. It has 
its source at the gaussian fixed point and stops at the 
Wilson-Fisher fixed point (it takes an infinite RG time 
to reach this fixed point). Many small rivers, coming 
from the mountains, fiow very rapidly "into" the large 
one. Each of them corresponds to a different initial con- 
dition of the RG flow, that is to a different microscopic 
system at scale More precisely, different "natural" 
initial conditions of the RG flow at scale A evolve at scale 
A: ^ A towards the same RG trajectory — the large river 
— up to an error of order j Is? which is the thickness, 
around the large river, of the beam of RG trajectories 
emanating from the set of natural initial conditions un- 
der consideration. 

Let us now consider a RG trajectory T emanating from 
a given initial condition at scale A. If at scale A: <C A, the 
difference between T and the large river L is neglected, 
the RG flow on T seems to be driven by a unique cou- 
pling that therefore determines, in the massless case, the 
whole theory. Within the perturbative scheme, the infi- 
nite cut-off limit enables to get rid of any reference to the 
cut-off scale, that is to A, and leads to a unique RG flow 
which is valid at any scale, at least in principle. It must 
therefore be the flow on L and the infinite cut-off limit 
has thus removed the transient regime where T and L 
are far away. The theory has thus been "projected" onto 

^^One should remember that the space of coupling constants is 
infinite dimensional and that it is therefore non-trivial that the 
projection of L onto the g4-axis be non vanishing. A randomly 
chosen vector in an infinite dimensional space has in general a 
vanishing projection onto a given direction. 

^^For the Ising model, different natural initial conditions could 
correspond for instance to different kinds of lattices or to different 
types of couplings among the spins: next nearest neighbor cou- 
plings, anisotropic couplings, etc. 



the line L (in the massless case). This is of course tech- 
nically convenient since the transient regime depends on 
the initial condition of the RG flow at scale A and thus on 
all couplings either renormalizable or nonrenormalizable 
involved in the initial condition. However, by taking the 
infinite cut-off limit, that is by neglecting the difference 
between T and L — which is a very good approximation 
in the IR regime — , one has given up the possibility of 
reversing the flow in the UV direction to go back on T 
to the initial condition at scale A. This means that it 
is no longer possible to initialize the flow at A from a 
"microscopic" action once T has been replaced by L (in 
the IR). To be consistent, the RG flow must therefore be 
initialized at a flnite scale, say /Lt, which is infrared-like 
with respect to A0 This is the role of the "renormaliza- 
tion prescriptions" of perturbation theory that enable to 
parametrize the theory in terms of a renormalized cou- 
pling (and mass) defined as the value of the renormalized 
correlation function at scale /i: 

rk'^({pf})=54,fi(M) (2.77) 

with {pf^} that are functions of /i. For instance, at a 
symmetric point 

.pf = /.2(%-l/4). (2.78) 

In the perturbative scheme, the renormalization pre- 
scriptions come from the necessity of fixing the arbitrari- 
ness of the finite parts of the subtraction procedure of 

Let us mention here that within dimensional regularization, 
for instance, there is no explicit ultraviolet regularization scale. 
The only scale introduced in this regularization scheme is the scale, 
often called /i, necessary to preserve dimensional analysis when, 
for loop-integrals, / d^g is replaced by J d'^q: / d^g — > /i*""* J d^q. 
This scale is also used most of the time as the scale of the renor- 
malization prescriptions that are either explicit, see Eg. 12.7711 . or 
implicit as in the MS scheme (the fact that they can be implicit 
does not change anything to our discussion). However, any regu- 
larization consists in modifying the short distance behavior of the 
theory because this is where the divergences come from. Thus the 
ultraviolet cut-off, although not explicit must be built from and 
e = 4 — d. Prom a comparison of the divergences obtained in di- 
mensional regularization and in the cut-off regularization, it is easy 
to get a qualitative correspondence between the two schemes and 
thus an estimate of the UV cut-off scale of dimensional regulariza- 
tion. When a logarithmic divergence occurs in the cut-off regular- 
ization scheme, a pole in e occurs in dimensional regularization. 
Thus it is reasonable to imagine that log A ^ To make this 

correspondence dimensionally valid we must have logA//i ~ 1/e. 
We therefore find that in dimensional regularization, the UV scale 
behaves as A ~ /iexp(l/e). 
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the divergences. Since oo + anything finite = oo, sub- 
tracting a divergence of a perturbation expansion is a 
well-defined procedure up to a finite part. To finally 
get a unique and well-defined renormalized theory, it is 
necessary to specify this finite part. Technically, this 
is what renormalization prescriptions do. From a RG 
point of view, they also enable to initialize the RG flow 
at scale Of course, and this can be checked on each 
example, there are as many independent parameters in 
a field theory as primitive divergences — and thus as 
renormalization prescriptions, let apart the field renor- 
malization — and this is also the number of dimensions 
of the attractive submanifold of the Wilsonian RG flow. 

Let us finally mention that the back-reaction of the 
nonrenormalizable couplings on the fiow of the renor- 
malizable ones can be analyzed within perturbation the- 
ory. In dimension four for instance, a gQcffi-teim in the 
action contributes to the fiows of the mass and of 54 in 
the following way. Since a new vertex exists, new graphs 
appear in the loop-expansion of r*^^-* and T^'^\ From 
the power counting point of view, that is by dimensional 
analysis, it is clear that these new graphs contribute to 
r^^) and r^'') as the (/)'*-coupling does, that is by factors 
and log A for F^^) and by factors log A to F^"*). This 
can be checked on the following graph contributing to 

r(2): 

Let us notice that the "divergences" coming from the 
integrals are more severe for the graphs involving a 
(ffi-vertex than for the others but that the overall A- 
dependence is indeed what is expected since ~ 
from dimensional analysis. For graph 12.151 for instance 



Initializing this flow in the IR has several advantages. First, 
as we already mentioned, this is anyway obligatory once the infi- 
nite cut-off limit has been taken since, in this case, any reference to 
a microscopic model defined at an UV scale has been lost. Second, 
if the model under study is not derived from a more fundamental 
model at scale A — in which case the analytical form of Fa is 
not known — , its initialization at scale A would require in general 
infinitely many phenomenological input parameters since Fa is in 
general not polynomial. This is of course impossible and should 
be compared to what is done in the IR: only the values of the 
renormalizable couplings have to be fixed since the dimension of 
the submanifold L is the number of renormalizable terms. For the 
Ising model for instance, it is possible to initialize the flow at scale 
A since already at this scale it is an effective model derived from a 
more fundamental model (coming for instance from the Hubbard 
model). The same occurs for fluids for instance: microscopic mod- 
els can be derived from other models that are more fundamental. 



Figure 2.15: The double tadpole diagram contributing 
to F*^^^ and that involves at first order the ge coupling. 
Each tadpole contributes to a factor A^ and the vertex 
(that is ge) to A^^. 

the two tadpoles contribute by a factor A"' and the ver- 
tex by a factor A~^. These new divergences, being of the 
same type than the ones coming from the 0*-coupling, 
can be eliminated through a renormalization of the field, 
the mass and of g^. The fiow of the renormalized 0'*- 
coupling, (?4,i{, is not affected by these contributions. 
From the RG point of view, this comes simply from the 
fact that once T and L are very close, the flow on both 
trajectories are almost identical: the difference between 
both flows that comes from a difference in the nonrenor- 
malizable couplings at scale A, shows up only in the first 
part of the fiow when T and L are far away from each 
other. More precisely, let us consider two initial condi- 
tions at the same scale A differing only by nonrenormal- 
izable terms, for instance two points Mi_a and M2.A on 
T and L corresponding to two different values of 55, A ■ 
At scale fc <C A these points have evolved in Mi^k and 
M2,k, see Fig |2.16l Although Mi^k is close to the trajec- 
tory L, the two points are in general not close to each 
other since the beginning of the two fiows are very differ- 
ent. This is why non universal quantities are in general 
not correctly computed if nonrenormalizable terms are 
neglected. However, one can consider the point M2 on 
L which is the closest to Mi^k and ^ the ancestor at 
scale A of this point. The two models described respec- 
tively by Mi^A and ^ differ at short distance but are 
(almost) identical at all scales larger than the scale fco 
where T and L are very close since all couplings in both 
models are almost identical for k < kg. In particular, all 
quantities — universal or non-universal — that can be 
computed from correlation functions at vanishing exter- 
nal momenta are found identical in both models. This 
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Figure 2.16: Flows on the critical surface for two differ- 
ent initial conditions Mi^a and M2.A € L differing only 
by the initial value of an irrelevant coupling. During the 
transient regime on T, that is before T is close to L, the 
RG flows differ between the two trajectories so that at 
scale k the points Mi ^ and M2^k are not close to each 
other although fc <C A and Afi ^ is almost on L. 

shows that, at least in principle and for sufficiently long- 
distance physics, both universal and non-universal quan- 
tities could be computed from the flow on L — and there- 
fore perturbatively — provided that one knows where to 
initialize the flow. Of course, the problem is that in 
general we do not have such an information and this is 
the very reason why non-universal quantities cannot be 
accurately computed within perturbation theory. 

2.6.2 Universality 

As explained previously, all the universal features of a 
given microscopic model can be computed from the flow 
on the large river. However, once the true RG trajectory 
T of the model has been approximated by the large river 
L in the long distance regime it is impossible to reverse 
the flow in the short distance direction to go back to 
the microscopic model we started with. Thus, once close 
to the large river, it is almost impossible to know from 
which small river comes the water, that is from which 



initial condition does the flow come from (this would re- 
quire high precision measurements). This is universality. 
We thus observe that universality is a much more general 
concept than the statement that some critical exponents 
are independent of the microscopic details. Critical ex- 
ponents are properties attached to the flow around the 
infrared attractive fixed point, if it exists. The attractive 
character of the submanifold that can be parametrized 
by renormalized couplings only is independent of the ex- 
istence of such a fixed point. It exists for theories that 
are not critical and it takes place far before the theory 
is in the deep infrared regime. It is this property that 
makes field theory — or at least perturbative field the- 
ory — interesting since this is what allows us to focus 
on a small number of couplings. The price to pay is that 
nonuniversal features are not (accurately) computable 
within perturbative field theory. 

Let us emphasize at this stage that the NPRG, because 
it is functional in essence, can follow any RG trajectory 
in the whole space of coupling constants contrarily to 
perturbation theory. The infinite cut-off limit A — > 00 
performed in perturbation theory and which is the short- 
cut that enables to eliminate any transient regime and 
to pick up the flow on the large river is not necessary in 
the NPRG framework. It is therefore possible to start 
at scale A with a given initial condition — polynomial 
or not in the field — and to integrate the RG equa- 
tion down to A: = 0. In principle, the full free energy 
of the model can be obtained in this limit and the only 
limitation of the method comes from the truncation of 
r[M] that has been chosen to integrate the RG equa- 
tion. Thus, contrarily to perturbation theory, there is 
no intrinsic limitation of the NPRG approach as for the 
computation of nonuniversal quantities. Several compu- 
tations of nonuniversal quantities such as critical tem- 
peratures in fluids [59] or phase diagrams for reaction- 
diffusion systems (the so-called branching and annihilat- 
ing random walks) |60j have been obtained accurately by 
NPRG methods. 

2.6.3 Continuum limit and asymptotic 
freedom 

Taking the continuum limit of a field theory consists in 
removing any reference to an ultra-violet cut-off while 
maintaining the long-distance physics fixed. At first 
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sight, this question seems equivalent to that of the infi- 
nite cut-ofF hmit addressed in the perturbative context. 
The difference conies from the fact that the question of 
the existence of the continuum hmit is addressed beyond 
perturbation theory. We shall only study the critical the- 
ory in the following although the massive case could be 
studied along the same line. Wilson's RG is the right 
framework to study this question. 

We shall study in the following the particular case 
where the "continuum limit", A — * oo, is taken on a 
particular RG trajectory, that is for a given physical 
system. For the massless Z2-invariant model in d = 3, 
there is only one trajectory for which this limit exists: 
the large river joining the gaussian to the Wilson-Fisher 
fixed point. The limit on this trajectory is taken in the 
following way. We imagine that an external information 
is known about the system at a scale k (an experiment 
has been performed for instance) that allows us to de- 
termine the value, at this scale, of at least one coupling 
constant. We are thus able to associate a point on the 
line L with the scale k. Keeping fixed this point, we can 
determine which ancestor at scale A corresponds to this 
point by inverting the RG flow in the ultraviolet direc- 
tion. Since we have chosen the large river as RG trajec- 
tory, the flow slows down as A is increased since the point 
associated with this scale {M2,a in Fig l2.16p gets closer 
and closer to the gaussian fixed point. When A — > oo the 
point associated with this scale tends to the gaussian and 
thus the continuum limit exists on this trajectory. The 
existence of a continuum limit is therefore a consequence 
of the fact that the theory is asymptotically free in the 
ultraviolet regime. One should be careful that, of course, 
this procedure should be taken as a limit: if one starts 
the RG evolution right on the gaussian fixed point, the 
representative point of the model remains gaussian. 

Let us emphasize that for any RG trajectory other 
than L, the continuum limit does not exist (in the mass- 
less case) since the flow blows up at infinity in the ultra- 
violet direction for these trajectories. Let us remind the 
reader that the flow is extremely rapid away from the 
large river. This means that as soon as the RG tra- 
jectory is not "close" to L, the flow in the ultraviolet 
direction diverges very rapidly from L. This does not 
mean that the corresponding fleld theory does not exist. 
This simply means that the ultra-violet cut-off cannot 
be removed since the flow is not controlled in the UV 



direction. However, as an effective (cut-off) field theory 
it is validin the IR and its universal behavior is the same 
as the field theory defined on L. 

Let us emphasize that the continuum limit is probably 
irrelevant from a physical viewpoint: who cares about 
the physics at asymptotically small distances or, in par- 
ticle physics, at asymptotically large energies ? This 
question goes beyond physics since, by definition, there 
will never exist any experiment able to test the infinites- 
imally short distance regime of any model. Field the- 
ory and more precisely renormalizable field theories are 
most probably only effective long-distance models for the 
underlying and more fundamental theories that proba- 
bly involve nonrenormalizable terms and/or are not even 
field theories. 

2.6.4 Perturbative (non-) renormaliz- 
ability again 

It is important at this stage to realize that the infinite 
cut-off limit that can be taken order by order of the per- 
turbation expansion of a renormalizable theory does not 
make reference to any RG trajectory and thus neither to 
the existence of fixed points nor to asymptotic freedom. 
It depends only on the possibility of removing recursively 
the divergences of the theory by the redefinition of the 
parameters of the model. We could thus question the 
meaning of this infinite cut-off limit. Let us first make a 
remark here. 

The existence of the infinite cut-off limit docs not im- 
ply that of the continuum limit. The example of the 
Z2-invariant theories in d = 4 is illuminating in this re- 
spect. When d is increased from 3 to 4, the Wilson- 
Fisher fixed point moves towards the gaussian and, in 
d = 4, coincides with it. The line L does no longer ex- 
ist in this dimension. The gaussian fixed point becomes 
infrared attractive and ultraviolet repulsive. Any point 
on the critical surface in d = 4 is attracted towards the 
gaussian in the infrared and pushed away from it in the 
ultraviolet. There is no continuum limit of any model 
apart from the gaussian model itselfF^ However, the 

^® Reciprocally, any model defined at an UV scale A on the crit- 
ical surface behaves at large distance as if it was free (this is the 
so-called "triviality problem" of the 0* theory in d = 4). This 
is the reason why d = 4 is the upper critical dimension of the 
Z2-invariant models. 
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infinite cut-ofF limit exists at all orders for the the- 
ory even in d = 4! This means that order by order, it 
is possible to construct finite renormalized Green func- 
tions defined by their series expansion in the (renormal- 
ized) (/)^ coupling without any reference to an UV cut-off. 
Of course, such an expansion does not guarantee that 
the sum of the perturbative series of the Green func- 
tions exist. Let us remind the reader that, as usual in 
physics, the perturbative series are not absolutely con- 
vergent. They are asymptotic series, at best. [ITl [CT [5] 
For the three-dimensional cj)'^ theory, it has been rigor- 
ously proven that these series are Borel-summable (in the 
"massive zero momentum scheme" their Borel-transforms 
are convergent). It is very probable that in d = 4 the 
renormalized perturbative series are not summablc for 
any finite value of (74, _r. Thus, the renormalized series of 
the Green functions can well exist at all orders, this does 
not mean that the functions themselves exist. 

It could seem strange to the reader that the pertur- 
bative RG is in a better position than the perturbation 
expansion of the Green functions to study the problem 
of the existence of the continuum limit. It is important 
to have in mind that, even at one-loop, the /3-function 
of the 0"* coupling constant: (3{gA,R) takes into account 
the leading-logarithms at all orders of the perturbative 
series. [6] The RG goes, in this respect, beyond perturba- 
tion expansion even if the P function is computed petur- 
batively since it is sensitive to informations on the be- 
havior of the sum of the series of the r(")'s that are 
unreachable at any finite order of their perturbation ex- 
pansion. This is the reason why, even at one-loop, the 
behavior of the (3 function of the 0'* theory leads to the 
conclusion (probably valid beyond perturbation theory) 
that the </)^ theory has no non-trivial continuum limit in 
d = 4 whereas the perturbation expansion of the Green 
functions exists at all orders at the infinite cut-off limit. 

Having in mind this subtlety, we could thus re-ask the 
question: what does mean perturbative renormalizabil- 
ity? We have already seen above that the infinite cut-off 
limit removes all finite cut-ofF effects from perturbation 
theory (the transient regime before reaching the large 
river L) and thus enables to pick up the RG flow on the 
line L. It thus allows us to select the (renormalizable) 
coupling(s) that drives the flow of all the others on this 
line. But then the relevant question is: what does mean 
perturbative non-renormalizability? After all, we have 



just learnt that all couplings are allowed (they are all non 
vanishing on the line L) and that they are all functions of 
the renormalizable one(s) once the RG trajectory has (al- 
most) collapsed with Thus, all couplings are present 
even in perturbation theory: they are simply not all inde- 
pendent. What does mean in this case that a coupling is 
perturbatively non-renormalizable? To answer this ques- 
tion we have to go back to what occurs in perturbation 
theory when a non-renormalizable coupling is present in 
the action of the model. In this case, at one-loop or- 
der, other non-renormalizable terms must be introduced 
in the action to cancel the divergences coming from this 
term. Then, at the next orders, these new terms gener- 
ate new diagrams and new divergences that themselves 
require new non-renormalizable terms in the action to 
cancel them. This process never stops and the total ac- 
tion involves infinitely many terms. This is not in itself 
sufficient to make the theory useless: we have already 
seen that in fact all couplings are allowed and are likely 
to be there at the scale of the cut-ofF A. The very point 
is that, as we already mentioned above, there must be as 
many renormalization prescriptions as there are indepen- 
dent divergences to fix the arbitrariness of the finite parts 
of the counter-terms. For a non-renormalizable theory, 
this means that there are infinitely many prescriptions 
and the theory is thus non predictive (each prescription 
requires an input, an experimental result for instance). 
This is the very reason why non renormalizable terms 
are supposed to be forbidden perturbatively. We can re- 
formulate and understand the non-renormalizability cri- 
terium from a RG point of view. Imposing a renormaliza- 
tion prescription on a Green function F*^") means impos- 
ing the value of this function for a given configuration of 
its momenta, that is specifying the value of the coupling 
gn at an infrared scale, say fj,. Let us now recall that 
the renormalization prescriptions are imposed indepen- 
dently of the RG flow within perturbation theory. The 
non-renormalizability of a theory simply means that it 
is "impossible" to impose that the RG flow of the theory 
goes at the infrared scale fi through a point arbitrarily 
chosen in the space of coupling constants: the flow must 
be on (or extremely close to) the large river L in the 
infrared. If we nevertheless impose that the flow goes 



^'^Once again we deal only with a massless theory. The presence 
of a mass does not change qualitatively the discussion below. 
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through such an arbitrarily chosen point, we thus find 
that (i) the theory is non-predictive even in the infrared 
since the flow is very rapid away from L and that there- 
fore a small change of scale below ^ — or a small error 
in the determination of the renormalization prescriptions 
— leads to very large changes of the running couplings, 
and (ii) the ancestor in the ultra-violet of the point cho- 
sen for the prescriptions is extremely far away from the 
gaussian: at A, the initial condition of the RG flow is 
thus extremely unnatural in this casein Thus, what 
perturbation theory calls non-renormalizable is a theory 
that is non-predictive in the infrared and non-natural 
in the ultraviolet. The subtlety here is that the reason 
why non-renormalizable theories should be discarded is 
not that it is forbidden to consider non-renormalizable 
couplings at scale A (they are in general present at this 
scale) but that it is impossible to choose at will their 
value in the infrared! The only possible choice to avoid 
the above-mentioned problem would be to pick up with 
the (infinite number of) renormalization prescriptions a 
point on L (or very close to L). In this case, the ancestor 
in the ultraviolet of this point would be also on L and we 
would have thus built at scale A a theory that pertur- 
batively would involve infinitely many couplings, would 
remain predictive in the infrared and well controlled in 
the UV since it would become gaussian asA ^ oo. In 
fact this theory would not really involve infinitely many 
couplings since as already emphasized, on L the flows of 
all couplings are driven by those of the renormalizable 
ones. This program has been initiated by Symanzik and, 
in Wilson's words, this is the best that we can do within 
perturbation theory. 



Let us now come back to physics and study the 0{N) 
models at order 0{d^) of the derivative expansion. 



■^*One should remember that we are deahng here with dimcn- 
sionless coupUng constants. In a natural theory, we expect all 
dimensionful coupling constants, at scale A for instance, to be of 
order of A. This precisely means that the dimensionless coupling 
constants must all be of order one. Of course, the continuum limit 
does not exist for a non-renormalizable model. 
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The 0{N) models at 0(92) of 
the derivative expansion 



Although the Ising and 0{N) models have much in com- 
mon, there are several non-trivial points specific to the 
0{N) models that are worth studying. Among them is 
the presence of Goldstone modes in the low temperature 
phase, the Mermin- Wagner theorem in d = 2 and the 
Kosterlitz-Thouless transition for = 2 in d = 2. We 
shall use again the derivative expansion that writes at 
0(92): 



(2.79) 



where M is a A^-component vector. In fact, we shall 
mainly study the LPA' that consists in neglecting 
Yk{M'^) and keeping only the first term of the field- 
expansion of Zk{M'^) that we call Zk (not to be confused 
with the partition function): 



1 



rx{Uk{M\x)) + -Zk{'^M) 



(2.80) 



The RG equation on Vk is obtained along the same 
line as in the Ising case. We start by constructing the 
partition function Zk[B] 



B. 



Zk[B] = J V(t>ix) exp i^-H[4,] - AHk[ 
with 



(2.81) 
(2.82) 



Since rfe[A/] is an 0(A^)-scalar, the RG equation involves 
now a trace on the 0{N) indices 



dk^k ^ Tr 



where (x, y) is now a N x N matrix: 

6^Wk 



dkRk{x-y)wlf\x,y) (2.83) 



Wt%{x,y) 



6B,{x) 5Bj{y) 



(2.84) 
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.(2) 



TZk is again the inverse of 



(2) 



5{x ~ z)5,u = j W^^lix, y) (4'^ + 7^fe)^^ (y, z) . 

(2.85) 

where repeated indices are summed over (Eintein's con- 
vention). Thus the RG equation writes: 



1 



(2) 



TZk] . (2.86) 

x,y 



2.7.1 The RG equation for the potential 

Once again we define the potential as rfe[M] evaluated 
in a uniform field configuration AI. By symmetry, we 
can choose any direction for M. We take 







M 



(2.87) 



V J 

The RG equation on the potential writes 



+ iZ,,q^ + Rk)S 




where the trace means summation over 0{N) indices and 
integration on q. Since 



■ + 



dMidMj dp ' dp- 
where p = 1/2M^, we obtain 



2 M^MJ 



(2.89) 



dM, dMi 



( Zkq^ +Rk + Ul + 2pU'^ 

Zkq^ +Rk + U'^ 



V 



2.90) 



Zkq^ + Rk + U'J 



It is trivial to invert this matrix and to compute the 
trace. We find: 

1 

2 



dkUk 



dkRkiq) 



Zkq^+Rk + U'^ + 2pU'^ 
N -I 



Zkq^ + Rk + U' 



(2.91) 



We can see two differences with what we have found in 
the Ising case within the LPA: 

• Within the LPA there is no "field renormalization" 
Z/j. in front of the q^ term. Its presence will have 
important consequences both from the technical and 
physical points of view. 

• There is a new term proportional to — 1 in 
Eq. (|2.9ip . It will take care of the physics of the 
Goldstone bosons in the low temperature phase. 

2.7.2 The RG equation for the dimen- 
sionless potential Uk 



Once again, Eq. (|2.9ip is not well suited for the search of 
fixed point since k appears explicitly in the right hand 
side through Rk and Zk- We have first to go to dimen- 
sionless quantities. But now, even the change of variables 
to dimensionless quantities, Ea. (|2.59p . is not sufficient to 
get rid of the fc-dependence since Zj. depends on fc. It 
can be shown that in the scaling regime ^52j 



(2.92) 



where 77 is the anomalous dimension. Thus, Zk never 
reaches a fixed point value and it is therefore necessary 
to get rid of it to find the fixed point potential. We there- 
fore introduce dimensionless and "renormalized" quanti- 
ties defined by 



y 



(2.93) 



Rk{q^) = Zk q^r (y) = Zk k^yr{y) (2.94) 
S; = kx (2.95) 
M(i) = ^/%k^ M{x) (2.96) 

Ut (M(i)) =k~'-Uk (M(x)) . (2.97) 

Note that a factor Zk has been included in Rk- We 
can now repeat all the different steps leading to the RG 
equation on Uk- It is useful to define first 



kdkZk 



-VkZk 



(2.98) 
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where r]k could be called a "running" anomalous dimen- 
sion. Because of the behavior of Zk , Eq. (|2.92p , rjk reaches 
a fixed point value (the anomalous dimension) whereas 
Zk does not. From Eqs.^^M^M^ we deduce that both 
Zkq"^ + i?fc + + 2pUl! and dtRk{q^) become propor- 
tional to Zkk^ . Thus, with this rescaling, the explicit k- 
and Zfe-dependences disappear in the equation for t/t-F^ 
The RG equation on Uk writes: ^32; 

dtUt = -dUt + {d-2 + rik)pU'^ 

poo 

-2vd dyy''/^ {mr{y)+2yr'{y)). 
Jo 

( I + ^ ^2 99) 

\y{l + r{y))+Ul + 2pU[' y{l + r{y)) + U[ ) 

We shall see in the following that it is sometimes con- 
venient to consider the field-expansion of Ut{p). In this 
case, the RG flow on the potential is projected onto an in- 
finite hierarchy of ordinary differential equations for the 
evolution for the coupling constants appearing in this 
expansion. It is a non-trivial question to know around 
which field configuration p one should perform the ex- 
pansion. Of course, if the field expansion of Ut{p) is 
not truncated and if its radius of convergence is infi- 
nite, the point around which the expansion is performed 
does not matter. However, the radius of convergence is 
in general not infinite and we shall be of course inter- 
ested in truncating the series expansion at orders in p as 
low as possible. The rule of thumb is that each time a 
field-expansion has to be performed, the best choice is to 
perform it around the minimum of the (dimensionless) 
potential Kk, Ea. (|2.65p : 



Ml 



(2.100) 



with 



V J 



(2.101) 



^®In this sense, working with the dimensionless and renormal- 
ized quantities consists in going to a "co-moving frame" where the 
exphcit fc-dependence has been eUminated. 



Let us notice that once Uk has been truncated at a finite 
order in p, this equation is not sufficient to define com- 
pletely Kk, XkjUs^k, etc. It is necessary to define them 



dUk 

dp lp=«fc 
d^Uk 



dp^ 



Qpn lp=« 



(2.102) 
(2.103) 

(2.104) 
(2.105) 



Note that Eq. p.l02p makes sense only if 7^ (for 
the search of the fixed point, this is not a problem since 

K* ^ 0). 

The flow of all these coupling constants can be ob- 
tained trivially by acting on both sides of these equations 
with dt and by using Eq. (|2.99p . We find for instance pS] 

dtKk = -{d-2 + r^k)Kk + 2vd (3 + 2^^^^ ) lf{2Kk\k) 



+2vdiN - l)lfiO) 



(2.106) 



dtXk = (rf - 4 + 2r?fe)Afc + 2vd{N - l)Xl4{Q) 



+2vd{3Xk + 2KkU3,kfl2{2KkXk) 



'2vd 2u3,k + 2KkUi.k 



l—^ lt{2KkXk) (2.107) 
At / 



where the so-called threshold functions /fj are defined in 
the Appendix, section [3.14l One remarks that the flow of 
Xk involves us^k and U4,fc. This is a general rule: the flow 
of u„^k involves Un+i,k and Un+2,k- The non-perturbative 
character of these flows comes from the non-polynomial 
character of the threshold functions l^. This, in turn, 
implies that the right hand side of Eqs. (|2. 10612. 107p are 
not series expansions in the coupling constant A^. 

The computation of the anomalous dimension rjk re- 
quires the computation of the flow of Zk- As we did 
for the potential this is possible only after a definition 
of Zk in terms of Tk has been found. It is clear that 
Zk corresponds to the term in Tk which is quadratic in 
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M and in q. In fact, this definition is not sufficient to 
completely characterize Zk since it is the first term in 
the expansion of the function Zk{p) and it is necessary 
to specify around which value of p the expansion is per- 
formed. Here again, we choose the minimum Kk of the 
potential, Eqs. ipTOTI - l^lU^ . 

A precise calculation shows that 

-(2) 

where ^\^2 p) (2 -p) second derivative of Tk with 

respect to M2{p) and M2(— p). The flow of Zk is now 
obtained by acting on both sides of (|2.108p with dt ■ After 
a straigthforward although somewhat tedious calculation 
we obtain: pS] 

Vk = ^^k \l ™^,2(2KfcAfc) (2.109) 

where 2 is a threshold function defined in (|3.17p . We 
thus find that rjk is not an independent quantity. It is 
entirely defined by the other couplings. 

2.7.3 The limits d ^ A, d ^ 2 and 

00 and the multicritical Ising fixed 
points 

Let us first recall that the whole perturbative series of the 
correlation functions can be reproduced from the exact 
RG equation on if it expanded perturbatively. Of 
course, our aim is not to use this approach perturbatively 
and it is thus interesting to understand the interplay 
between the derivative expansion and the perturbation 
expansion. 

A very nice feature of the effective average action for- 
malism is that the one-loop results obtained in d = 4 — e 
(with the 0^ theory) and d = 2 + e (with the non lin- 
ear sigma model for iV > 3) are retrieved very simply 
while keeping in the ansatz for Tk only Kk , Afc and Zk , 
Eas. l|2. 10612. 10712. 109p . [5 ^ [55] This means that the (3- 
functions cither of the coupling Afe around d = 4 or of 
the temperature (which is related to the inverse of Kk) 
around d = 2 reproduce at leading order the one-loop 
result respectively in d = 4 + e and in d = 2 + e. The 
large N limit is also retrieved at leading order with the 



same ansatelf^ This is a very interesting property since 
it shows that the same calculation leads to controlled 
results both in the upper and lower critical dimensions 
and at = oolf^ The NPRG results are thus, at least, 
a clever interpolation between the perturbative results 
obained in = 4, d = 2 and at A^ = 00 ! Let us 
emphasize that this should be enough to convince the 
most reluctant that the NPRG leads to results that are, 
at least, not less controlled than the perturbative ones. 
They are in fact probably much better controlled since 
they are one-loop exact in the two critical dimensions 
and since the three-dimensional critical exponents seem 
to converge rapidly with the order of the derivative ex- 
pansion, see section [2?T] 

Of course the A" = 1 and N = 2 cases are peculiar 
since for A^ = 1 the lower critical dimension is one and for 
N = 2 there exists a finite temperature phase transition 
in d — 2 of inifinite order: the Kosterlitz-Thouless tran- 
sition induced by the vortices. As for A^ = 1, the order 
0{d^) approximation leads to results in d = 2 that are 
qualitatively correct and quantitatively not so bad (er- 
ror on the critical exponents around 25%^). In d = 1 
the anomalous dimension is 1 and it is thus difficult to 
reproduce it within the derivative expansion!^ 

As for A^ = 2 in d 2, it has been shown — still 
with the ansatz involving only Kk , and Zk — that the 
Kosterlitz-Thouless transition is qualitatively well repro- 
duced. Stricto sensu, no line of fixed points is found but 
rather a line of quasi-fixed points where the flow does 
not stop but is very slow. On this "line" the correlation 
length is not infinite but very large. With the complete 
0{d^) approximation, a fairly good quantitative agree- 
ment is obained. [57] The remarkable point here is that 
these results have been obtained without introducing by 
hand the vortices as is usually done otherwise through 
the Villain's trick. This is very encouraging for the study 
of systems where the vortices play an important role but 

'"^Let us emphasize that this is not the case with the Wilson- 
Polchinski approach: even the one loop result in d = 4 — e is not 
reproduced at any finite order of the derivative expansion. [62] 

Needless to say that this is completely out of reach of the per- 
turbative expansions performed either from the 0'* theory (around 
d = 4) or from the non-linear sigma model (around d = 2). 

Although not proven, it is reasonable to assume that the 
smaller the anomalous dimension the better the convergence of the 
derivative expansion (»? = within the LPA). Thus r) is perhaps 
the "small parameter" of the derivative expansion. 
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where there is no spin waves/vortices decoupling and 
thus no simple analytical approach available. 

Morris has also shown that the infinite sequence of 
multicritical fixed points of the Ising model in d = 2 can 
be retrieved using the NPRG. This is also a non-trivial 
result since it would be very complicated to obtain them 
perturbatively. [56j 

2.8 Other fields of application of 
the NPRG in statistical me- 
chanics 

Apart from the 0{N) model, the NPRG has been used 
in the study of many systems both in particle physics 
|63j and in statistical mechanics. Let us mention a few 
of them in this latter field. 

• The magnetic frustrated systems. For some antifer- 
romagnetic systems the spins in the ground state 
are not coUinear. The system is then said to be 
frustrated. For instance, continuous spins on a tri- 
angular lattice adopt a 120° structure at T = 
when they interact antiferromagnetically since the 
spins on a triangular plaquette cannot be all anti- 
aligned. For three-component spins this 120° struc- 
ture implies that the rotational SO {3) symmetry is 
completely broken down in the ground state con- 
trary to the ferromagnetic case. Both the number 
of Goldstone modes in the low temperature phase 
and the critical physics are thus different for these 
systems and for non frustrated ones. The frustrated 
systems belong to the class of systems for which the 
symmetry breaking scheme is 0{N) 0{N — 2). It 
is amazing that their critical physics in d = 3 is not 
yet fully elucidated: the NPRG approach predicts 
that for N ~ 2 and TV = 3 in d = 3 generically weak 
first order transitions will be observed[35 whereas 
a fixed point is found perturbatively at five and six 
loops which implies a second order behavior [641 165] . 
Neither the experiments nor the numerical simula- 
tions enable to discriminate up to now between the 
two scenarios. 

• Out of equilibrium phase transitions. For systems 
that do not satisfy detailed balance the probability 



distribution of equilibrium states is not known a pri- 
ori: there is no analogue of the Boltzmann weights. 
It is however a matter of fact that many characteris- 
tic features of systems at thermal equilibrium exist 
also out of equilibrium. For instance many such sys- 
tems undergo continuous phase transitions governed 
by power law behaviors exhibiting universality. The 
whole machinery of field theory and renormalization 
group can be used with the subtlety that the "static" 
properties can be computed only through the large 
time limit of the underlying dynamics. For many 
systems this dynamics is given at a mesoscopic scale 
by a Langevin equation, the formal solution of which 
is given through a field theory. NPRG techniques 
have been adapted to this type of models. [66] Let 
us mention two non trivial results found this way in 
the so-called reaction-diffusion systems. 

One of the simplest system studied consists of par- 
ticles A on a lattice that can diffuse (with a rate 
D), spontaneously decay (with a rate n), annihilate 
by pairs when they meet on the same site (with a 
rate A) and give birth spontaneously to a daugther 
particle (with a rate a): A 0,2A ~* 0, A ^ 2A. 
The physics of the system consists in the compe- 
tition between the creation and the annihilation of 
particles. Depending on the magnitudes of /i, A and 
a the system either goes at large time into an ab- 
sorbing state where all particles have disappeared or 
to an active state where the average density reaches 
a non vanishing asymptotic value. In between these 
two situations exists a continuous phase transition, 
the universality class of which is famous: it is the so- 
called directed percolation universality class. How- 
ever, for /i = (no spontaneous decay) and ct ^ 0, 
the system is predicted (i) at the mean field level 
to always reach the active phase, (ii) from pertur- 
bative RG to undergo a continuous phase transi- 
tion between an absorbing and an active phase in 
d = 1 and d = 2 and only in these dimensions and 
(iii) from the NPRG to undergo a continuous phase 
transition in all dimensions. Numerical simulations 
have shown that a phase transition indeed exists in 
all dimensions and that the phase diagram found 
thanks to the NPRG is in quantitative agreement 
with the numerical data. [60j This is of course possi- 
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ble only because the NPRG is able to compute non 
universal quantities such as a phase diagram. For an 
other kind of reaction-diffusion system for which the 
parity of the number of particles is conserved by the 
dynamics {2A ^ 0, ^4 ^ 3A) the universality class 
is different and the fixed point was not satisfacto- 
rily found within perturbation theory. Thanks to 
the NPRG this has been achieved and the determi- 
nation of the critical exponents is in good agreement 
with the numerical data. [57] 

• Disordered systems. Disorder can be relevant for 
the critical physics of a system. Two kinds of dis- 
order have been studied for the Ising model: dis- 
order due either to randomness in the strength of 
the magnetic couplings Jy between spins (random 
mass Ising model) or to the coupling to an external 
random magnetic field (random field Ising model). 
The first type of disorder can be studied rather eas- 
ily by perturbative means and by NPRG methods. 
The net result of this type of disorder in three di- 
mensions is to sligthly modify the value of the crit- 
ical exponents. [68j The other type of disorder has 
a much drastic effect and is far subtler to study. 
Perturbatively, it has been proven at all orders that 
the critical physics of the disordered system in di- 
mension d is identical to that of the pure system in 
dimension d — 2 (this is the so-called dimensional 
reduction). Although true at all orders of pertur- 
bation theory, this result is wrong. The recourse to 
functional methods is unavoidable in this case since 
it can be shown that the effective potential devel- 
ops non-analyticities (a cusp) along the RG flow at 
a scale ko called the (inverse) Larkin's length. It is 
therefore necessary to renormalize this function it- 
self and not only the coupling constants that are its 
Taylor coefficients since the Taylor expansion ceases 
to exist below a scale fco of the RG flow. The func- 
tional RG as well as its nonperturbative version, 
the NPRG, has been successfully employed in this 
context. [691 EQl EH Ea [13 El] 

• Lifshitz point, bubble nucleation, Bose-Einstein con- 
densation. 

Let us finally mention other results obtained in sta- 
tistical mechanics with the NPRG. 



The study of Lifshitz critical points that is of critical 
systems exhibiting modulated phases (in space) has 
been achieved by Bervillier. 36j This requires to take 
into account anisotropic derivative terms as well as 
terms of order four in the derivatives. This kind 
of critical phenomena exhibits also genuinely non 
perturbative phenomena. 

Bubble nucleation is one of the most important phe- 
nomena occuring at a first order phase transition. 
The calculation of the nucleation rate is far from 
trivial and has been computed using the LPA. It 
compares well with other approaches. [75] [76] 

Recently, the problem of the cross-over between 
ECS superfluidity and Bose-Einstein condensation 
in fermionic systems has attracted much attention in 
particular because of experimental breakthroughs. 
The NPRG is ideally suited to study the cross-over 
between these phenomena since non universal quan- 
tities can be computed out of it and since it is 
functional. A quantitative comparison between ex- 
perimental data and theoretical predictions is now 
possible. [Ill EH] 

The NPRG has also been applied to the study of 
simple fluids. A study on the special case of CO2 
has been performed [59] with good results for the 
calculation of the critical temperature for instance 
and a general formalism has been developped in [791 
[80]. 



In this introduction to the NPRG we have focused on its 
application to statistical mechanics and on some of its 
relations with perturbative renormalizationlf^ We have 
seen three important points. 

First, at a conceptual level, the NPRG enables to un- 
derstand how microphysics can be continuously related 
to macrophysics, something that is not possible in gen- 
eral within perturbative field theory. As a by-product, 
one can solve this way the paradox that a field theory 
involving infinitely many interacting degrees of feedom 

^^For a pedagogical inroduction to the NPRG applied in gauge 
theories see ref. i63i and for a discussion of some features of the 
NPRG see ref. [51]. 



2.9 Conclusion 
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can be described in the infrared regime with only a fi- 
nite (and small) number of coupling constants, precisely 
those that are called renormalizable within perturbation 
theory. This comes from the attractive character of the 
submanifold spanned by the renormalizable couplings in 
the space of coupling constants (the large river effect) 
and is the very meaning of universality. 

Second, we have seen that contrary to common be- 
lief, it is possible to obtain qualitatively good results 
about the long distance physics with NPRG techniques 
from very short ansatze and even results that reproduce 
one-loop results around the upper and the lower critical 
dimensions and at large N. This is probably why the 
NPRG results obtained at finite N and for dimensions 
in between the upper and the lower critical dimension 
are reliable. 

Third, the series obtained from the derivative expan- 
sion seem to converge rapidly, at least in dimension three 
for the Ising model. This makes the NPRG a quantita- 
tive tool for studying strongly correlated systems and 
not only, as often claimed, a qualitative one. Of course, 
this claim shoud be substantianted by calculations per- 
formed beyond the 0(9^) and also in dimension two. It 
is nevertheless encouraging to see that critical exponents 
already converge at this order to the best known values 
without any resummation and that non universal quan- 
tities can be accurately computed. 

Let us finally mention that a crucial drawback of the 
derivative expansion is its inadequacy to the calculation 
of the momentum dependence of the correlation func- 
tions. In fact, it can be shown that the derivative ex- 
pansion makes sense only when the external momenta 
of the correlation functions are less than the running 
scale k. Thus, when k ^ only the infrared physics 
can be computed with the derivative expansion. Cru- 
cial improvements in the computation of the momentum 
dependence of F^^^ and F^^^ has been perfomed these 
last years [52l [53l [54] and there is no doubt that if this 
method works it will be a new step in our possibility 
of computing new non perturbative phenomena in field 
theory. 
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Chapter 3 

Appendix 



"Nobody ever promised you a rose garden." 

J. Polchinski 



3.1 Definitions, conventions 

• Integrals in x and q spaces 

In real and Fourier spaces we define 



• Correlation functions 

T[M] is a functional of AI{x). We define the IPI cor- 
relation functions by 



r(")[M(a;);a;i,...,x„] = 



SM^^ . . . SM^ 



Fourier transform 



(2^ 



(3.1) 



/ /(g)e'*- , f{q)^ f{x)e-^'>- . (3.2) 

J q J X 

Depending on the context, we omit or not the tilde on 
the Fourier transform. 



We also define the Fourier transform of F^"^ by 
f(")[Af(x);<zi,...,<7„] = 

4...x„e-^^'«'^'r(")[M(a;);xi,...,x„] . 
and also 

(5"F 



f(")[M(a:); 



, 9«] = 



(3.6) 



(3.7) 



(3.8) 



Definition of Vd 

f{q^)^2val dxx'^'-'fix) 



The relation between f ^"^ and F'"^ follows from Eq. 

f (") [M(x); <zi, . . . , g„] = (27r)-"'' f(") [M{x);q^, . . . , g„] . 

(3.9) 



with 



1 



Vd 



EYinctional derivatives 



(3.3) 
(3.4) 



• Cut-off function in x and q spaces 

1 



AHk[4>] 



Rk{x-y) (f>y 



5cj){q) 5(j){x) (27r) 



1 

(3.10) 

One should be careful about the fact that Rk is some- 
times considered as a function of q and sometimes as a 
function of q^ . It can be convenient to define a cut-off 
function with two entries by 



TZk{x,y) = Rk{x - y) . 



(3.11) 



52 



Then 



• Derivative of li 



nk{q,q') = {27r fS'{q + q')Rk{q) 



(3.12) 



dJ'iiw,Tj)^-{n + d„,o)lt+iiw,v) (3.20) 



• fc-dependent anomalous dimension 

By definition: 

kdkZk = -rjkZk ■ 

• Threshold functions l"^ 

n + (5„ ,0 



• (?-cut-off 

A convenient cut-ofF function Rk tliat allows to com- 
pute analytically some threshold functions is 



(3.13) 



Rk{q)^Zk(e-q^)9[l 



With this cut-off we find 



21 

fc2 



where 



dy y 



d/2-l 



s{y) 



(3.21) 



(3.22) 



(y(l + r{y)) + w 



n+l 



(3.14) 



Threshold functions and n with the 6'-cut- 



ofT 



Rkiq ) - Zkq r{y) with 2/ - p (3.15) ^Vith the cut-off function, Ea. dX^ . the 1^ threshold 

and, by definition of s{y) 



kdkRkiq') - kdk[Zkq^r(^ 



functions can be computed analytically 

Vk \ 1 



Zkk^~7jkyr{y)-2y\'{y)) 
Zkk'^siy) . (3.16) 



it{w,v) = ^{n + Sn,o) (^1 

777-2,2 (?«) = 



d + 2j (1 + w)"+ 

1 



{1 + wy 



(3-23) 
(3.24) 



• Threshold functions 



3.2 Proof of Eq.(l23I]) 

We define: 



,/2. {d^pf{x,Q) 



J = I e 



with 



P"i(a;,0)P"2 (a;,M;) 

(3.17) , . , 

and we rewrite the exponent: 



P{x, w) = Zk X + Rk {x) + w 



(3.18) y2 ,2 



• Universal value of Z^j"(0,0) for n > 

For n > and independently of the choice of cut-off 
function Rk'- 

e(0, 0) = f r dy (-2) = 1 (3.19) 

2 Jo (l + r(?;)) 



2a 2/3 2 ^ /3 2/3 



1 7 



We now define 



a 

u = y X 

1 



(3.25) 



(3.27) 
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and change variables: (y, z) — > (u,x). The jacobian is 1 
and thus: 



J = / e" 

/ U.X 



l2Tra(3 



(3.28) 



(3.29) 



3.3 The exact RG equations 

For the sake of simplicity, we consider a scalar theory 
(e.g. Ising). We have by definition 



Zk[B] ^ j'D<li exp l^~H[^] - Ai/fe[(/)] + ^ ^0 
with /^Hk[4>] = ^J Rk{q) 



"9 •r-q 



Wk[B]^ log Zk[B] 



Tk[M] + Wk[B] - / BM - i /" i?fe,:._j; Af, 

-/a; ^ J x,y 



with, by definition of M{x) 
6Wk 



5B{x) 



M{x) = {<j,{x)) 



(3.30) 



When B{x) is taken fc-independent (as in Zk[B]) then 
M{x) computed from Wk is fc-dependent. Reciprocally, 
if M(x) is taken fixed (as in rfe[M]), then B{x) computed 
from Eq. ()3.34|) becomes fc-dependent. 

3.3.1 RG equation for Wk[B] 

5fee'^'= = - i y" ( ^ 0, dkRk [x - y) 4>y) . 

. exp ( - H[4>\ - i y" Rk{q)M-q + / 



dkRk{x - y) 



S S 

sKJb^ 



(3.31) 



We therefore obtain for Wk- 

^ M/ [Rl ^ f f> TJ f A fJ^Wk ^ SWk 5Wk 

OkWkiB] ^ - ^ dkRki.^y) [j^^ + ^ ^ 

(3.32) 

which is equivalent to the Polchinski equation. 



3.3.2 RG equation for Tk[M] 



We first derive the reciprocal relation of Eq. (|3.30p . The 
Legendre transform is symmetric with respect to the two 
functions that are transformed. Here the Legendre trans- 
form of Wk is Tk + 1/2 J RkMM. Thus 

Tfc + i y Rk {x - y) = B^ (3.33) 



and then 



(3.34) 



In the Polchinski equation (|3.32p . the fc-derivative is 
taken at fixed Bx- We must convert it to a derivative at 
fixed M: 



dkMx 



\b 



dMx 



(3.35) 



Acting on Eq. (|3.30[) with dk\g, we obtain: 

dkrk[M]\^ + dkWk[B]^^ = / BdkM\g 

~\ I dkRk,x-y Mx My - / Rk,x-vMAMy^^M) 

^ J x.v J x.y 



Subsituting Eqs. (|3.34|3.32|3.35| into this equation we 
finally obtain 

dkTk[M]^\j^ dkRk{x-y)^^^ (3.37) 

The last step consists in rewritting the right hand side of 
this equation in terms of Tk only. We start from p.30p 
and act on it with S/SMz'- 

. S^Wk fPWkSBy 
S{x-z) ^ ^^^^ ^ , (3.38) 



SB, 6A'h Jy SBx SBy SM, 



54 



Now, using (|3.34p . we obtain the local potential approximation (LPA) that consists in 

truncating Tk as in (|3.44p with Zk{M) = 1: 

We define 



SB^ SBy \SMySM, ' -7 ■ rl^^[Mix)]= I i Uk{M^ (x)) + - (V Alf ] . (3.46) 



(3.39) 



2 



j2^^ By acting on Ea. l|3.45p with dk wc obtain: 



6B, SBy 



and thus '^^^ Jq ^ ' 



dMM) = ^J dMq){Tf^ u,_, 

(3.47) 

,5(x - z) = / T^i^^(a;,y) (r^^^ +7^fe) (y, z) . (3.41) Thus, we have to invert rf^ + Rk for a uniform field 
Jv ^ ^ confiffuration and within the LPA. From now on. we omit 



configuration and within the LPA. From now on, we omit 
the superscript LPA on F^. An elementary calculation 

F^^'' + TZk is therefore the inverse of wjf'' in the operator leads to 

(2) 

sense. Note that although we did not specify it, Wj, ' is f d'^Uk \ 

a functional of B{x) and T^^^ a functional of M{x). Re- '^kl,q' = \ -q]^ + j {"^^y^^il + ■ (3-48) 
lation (|3.4ip is valid for arbitrary M . The RG equation 

P-37P can now be written in terms of F^ only: Using 5{q = 0) = rj(27r)"'' wc find 

dkTk[M]^\( dkRk{x^y)(T'-^Unky\x,y) . dkUk ^ U ^^^^^^^^K^ ■ (3.49) 

(3.42) +^^('^)+aM^ 

In Fourier space this equation becomes: 

It is convenient to re-express this equation in terms of 



dkVkm = ]- I dkRkiq) (ff -f Ttfc) 



(3-43) _ 1,. 



P = (3.50) 



3.3.3 RG equation for the effective po- 



which is the Zo-invariant. 



tential 1 r dkRkjg) 

dkUkip) - / ^ \^TP( TT'i \ 3-51 

2 J + Rk{q) + Ul{p) + 2pU',:{p) 



The derivative expansion consists in expanding F^ as 

(3.44) 



-. ^ where Uk{p) and Uk{p) are derivatives of Uk with respect 



Tk[M{x)] = / ( C/fc(l'r) + -Zk{hn (VMY + ...] top 



To obtain the RG equation for the dimensionless po- 

, , , 1 , , , , 1 . . tential we have to perform the change of variables of 

where we have supposed that the theory is symmetric „ „„ , , 
, „ V f 4- f Ai^2 1 rr. Eq l2.60l We find 

so that Uk,Zk, ■ ■ ■ are functions of onfy. lo com- 
pute the flow of these functions it is necessary to define Q 
them from F^. The effective potential Uk coincides with ^^Mp = ^^^Ip + (d- 2)p ^ (3.52) 

Tk when it is evaluated for uniform field configurations 
Munif.: and 

rfe[M,„i,] = n Uk{Ml^si) (3.45) kdk\^. = kdk\^ " 2fc' ^'(v) (3.53) 

where is the volume of the system. It is easy to de- Inserting these relations together with F^q.^M in 
rive an RG equation from this definition of Uk if we use Eq-dSH]) leads to the RG equation on Ut, Eq.jMB). 
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